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ABSTRACT 



Context. Observations with space-borne X-ray telescopes revealed the existence of soft, diffuse X-ray emission from the inner regions 
of planetary nebulae. Although the existing images support the idea that this emission arises from the hot shocked central-star wind 
which fills the inner cavity of a planetary nebula, existing models have difficulties to explain the observations consistently. 
Aims. We investigate how the inclusion of thermal conduction changes the physical parameters of the hot shocked wind gas and the 
CIh' amount of X-ray emission predicted by time-dependent hydrodynamical models of planetary nebulae with central stars of normal, 

D ' hydrogen-rich surface composition. 

^0 I Methods. We upgraded our ID hydrodynamics code NEBEL by to account for energy transfer due to heat conduction, which is of 

importance at the interface separating the hot shocked wind gas ('hot bubble') from the much cooler nebular material. With this new 
version of NEBEL we recomputed a selection of our already existing hydrodynamical sequences and obtained synthetic X-ray spectra 
for representative models along the evolutionary tracks by means of the freely available CHIANTI package. 
ri [ Results. Heat conduction leads to lower temperatures and higher densities within a bubble and brings the physical properties of 

O 1. the X-ray emitting domain into close agreement with the values derived from observations. The amount of X-rays emitted during 

I ' the course of evolution depends on the energy dumped into the bubble by the fast stellar wind, on the efficiency of 'evaporating' 

O ' cool nebular gas via heat conduction, and on the bubble's expansion rate. We find from our models that the X-ray luminosity of a 

planetary nebula increases during its evolution across the HR diagram until stellar luminosity and wind power decline. Depending 
C/2 , on the central-star mass and the evolutionary phase, our models predict X-ray [0.45-2.5 keV] luminosities between 10"** and 10"* 

^^ ■ of the stellar bolometric luminosities, in good agreement with the observations. Less than 1% of the wind power is radiated away in 

this X-ray band. Although temperature, density, and also the mass of the hot bubble is significantly altered by heat conduction, the 
dynamics of the whole system remains practically the same. 
""►-C ■ Conclusions. Heat conduction allows the construction of nebular models which predict the correct amount of X-ray emission and 

y-^ • at the same time are fully consistent with the observed mass-loss rate and wind speed. Thermal conduction must be considered as a 

-«^ ' viable physical process for explaining the diffuse X-ray emission from planetary nebulae with closed inner cavities. Magnetic fields 

-vj I must then be absent or extremely weak. 

CO ' Key words, heat conduction - hydrodynamics - planetary nebulae: general - planetary nebulae: individual (NGC 2392, NGC 3242, 

h^ ! NGC 6543, NGC 7009, NGC 7027) - radiative transfer - X-rays: stars 
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^^ ■ 1. Introduction layers of t he PN. For recent reviews see [Schonberner & StefFenl 
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( l2003h and lSteffen & Schonbernei] (|200€ 
The modern, very successful concept for the formadon and evo- ^^ ^^^^^^^ ^.^^ ^^^ including the thin transition regime 



lutionofthe main structures ofPlanetary Nebulae (PNe) IS based ^^^^^^ ^j^^ p^ ^^^ ^j^^ j^j^^j ^-^^^ ^^ ^^^^^ X-ray 
H on the dynamical effects caused by (i) the interaction of a rapidly emission, mainly by thermal bremsstrahlung and line emission 
varying central-star stellar wind with the slow AGB wind ejected jVoik& KwoklllQSl . Any positive detection of diffuse X-ray 
eai-her, and (ii) the heating of this circumstellar material by pho- ^^^..^^^ coming from the inner cavities enclosed by nebular 
toionization. The wind from the ceiitral star is very fast exceed- ^j^^jj^ .^ ^ ^^^^^ confirmation of the wind interaction scenario 
ing in most cases 1000 km s"' and passes through a shock be- ^^ ^^^^^^^^ ^^^^^ ^^^^ ^^^^.^ ^^^ considerable interest to ob- 
fore making contact with the dense slow AGB wind A tenu- ^^^^^ ^^j^^^^^ p^^ ^-^^ X-ray satelhtes once they became avail- 
ous but very hot (^ 10^ K) bubble is formed since 9/16 of the ^^j^ ^.^^^ ^^^ surface bright ness in X-rays is expected to be 
kinetic wmd energy is converted into internal energy behind a .^^ ^^^jj ^see lVolk & KwoklfToSSh . only very few positive de- 
strong adiabatic shock. The bubble is separated from the outer, ^^^^-^^^ ^^^^ ^^^^ ^ ^^^^^ ^^ ^^^^ Summaries of X-ray ob- 
much cooler nebular gas by a contact discontmuity (or contact .e^vations of P Ne conduct ed to d ate with th e different satellites 
surface). The bubble s tiiermal expansion accelerates the mner ^^^^,^^ in lchu et al.l (l2()03h . lGuerrerol(l20()6l) . and lic^^li^ 

(I2007h . 



Send offprint requests to: M. Steff'en, e-mail: mstef f en@aip . de Although the observational confirmation of the existence of 

* Dedicated to the memory of M. Perinotto, a dear friend and es- shocked wind material was very gratifying, the properties of 

teemed colleague who died unexpectedly and much too early on August the X-ray emitting gas were disturbing. It turned out that the 

15, 2007. X-ray spectra are rather soft, indicating temperatures of only 



* (1 ... 3) X 10^ K. The electron densities of the emitting vol- 
umes vary between about 20 and 200 cm"-'. These results are in 
sharp contrast with theoretical expectations since the observed 
wind velocities and mass-loss rates demand bubble temperatures 
of X! 10^ ... 10^ K, and electron densities far below the observed 
ones. 



lAkashi et all (|2006) and lAkashi et al.1 ( |2007|) proposed that 
the remains of the much slower wind blown during the early PN 
evolution ('early wind') when heat conduction is unimportant 
are responsible for the properties of the observed X-ray emis- 
sion. Also collimated outflows (jets) can play a role. So it ap- 
pears that several physical mechanisms exist (heat conduction 
and/or mixing, early wind, jets) whi ch all may c ontribute to 
the X-ray emission, as pointed out bv lSoker & Kas tner (2003). 
The jet- wind interaction has recently been studied in detail by 
lAkashi et al.1 (I2008.) . 

In two fundamental papers on t he properties of inte rstellar 
bubbles by 'Castor et all (|1975') and ' Weaver et all (|1977[) it was 
demonstrated that heat conduction by electrons across the con- 
tact surface which separates the hot shocked wind gas from the 
much cooler swept-up matter is a natural mechanism to account 
for the observed X-ray and EUV line emissions. Heat conduction 
enforces 'evaporation' of cool gas into the hot bubble, leading to 
a shell of matter at the bubble's outer edge with properties just 
ideal for explaining the observed X-ray and EUV line emission. 

Since planetary nebulae are virtually scaled down versions of 
the interstellar bubbles (or H ii regions), heat condu ction should 
be important for them as well. To our knowledge. Soke? (1 19941) 
was the first to look into the consequences of heat conduction in 
PNe. He concluded that the observed X-ray emission must come 
from the heat conduction front, and not from the hot bubble as a 
whole. Heat conduction is also capable of explaining the unusu- 
ally strong UV lines (e.g. of O vi) seen in so me PNe whic h can- 
not be produced by regular photoionization. lSokeJ (1 19941) noted 
also that weak tangential magnetic fields would suppress heat 
conduction very effectively, and hence also the X-ray emission. 
But even without magnetic fields X-rays from old, large PNe 
may escape detection because the X-ray surface brightness falls 
below the detection limit of existing X-ray satellites. We note in 
passing that the same considerations apply to wind-blow n bub- 
bles around massive stars as well (cf. Wrigge e t alJll994l) . 

Despite of their potential importance, quantitative theoretical 
calculations of X-ray and EUV emiss ions from PNe based on the 
concept of heat conduction are rare. iZhekov & Perinottol (1 19961 
Il998h succ eeded in deve l oping analytical solutions based on the 
concept of Weave r et alj (11977') and were able to present quan- 
titative predictions for a fe w specified cases. Us ing analytical 
solutions of heat conduction. iGruendl et alJ (l2004l) modelled the 
transition layers between the hot bubble and the nebular regime 
for NGC 6543 in order to compute the O vi emission lines seen 
in the FUSE spectra. 

Two hydrodynamical studies addresse d recently the issue 
of X-ray emission f rom planetary nebulae dAkashi et al.l 120071 
IStute & Sahaill200 o). The basic philosophy of both studies is 
to avoid thermal conduction and to select the properties of the 
wind such that temperature and density of the X-ray emitting 
gas agree with the observations. No relation between stellar pa- 
rameters and wind properties is considered, and the feedback 
of radiation to the hydrodynamics is neglected. To explain the 
observed X-ray luminosities, these authors came up with wind 
parameters which disagree strongly with both the predictions of 
the theory of radiation-driven winds (e.g. Pauldrach et alJfl988l) 
and with the directly observed mass-loss rates and wind speeds 
of the objects in question. More precisely, the necessary mass- 
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loss rates are too large and the corresponding wind speeds too 
low - by as much as a factor SQ 

The results of l Akashi et al.l (l2007h are particularly interest- 
ing. Consistence with existing observations in terms of X-ray 
luminosity and temperature of the X-ray emitting gas could be 
achieved only for rapidly decaying central-star winds: while the 
wind speed increases linearly with time, the mass-loss rate must 
decline quite rapidly to very small values (see'Akashi et al .'20071 
Figs. 6 and 8 therein). For the best choice of Akashi et al. (2007|), 
the central-star wind reaches 10""^ Mq yr"' already at a modest 
wind speed of only 750 km s"' . 

We iterate that combination s of mass-loss rate s and out- 
flow speeds as they are found in lAkashi et al.l (12007!) are in se- 
vere conflict with both current theories of radiation-driven winds 
from hot stars (see Fig. |2] in Sect. [3] for details) and with the 
observations (see Table [3] in Sect. 15.21 ). with the apparent ex- 
ception of NGC 2392. In our opinion the existing studies of the 
X-ray emission from planetary nebulae which neglect thermal 
conduction by electrons are therefore far from being convincing. 
After all, thermal conduction is a physical process inherent to 
all hydrodynamical systems and becomes important in rarefied 
plasmas wherever the mean free path of the electrons is large 
enough. Thermal conduction can only be modi fied or even sup- 
ressed by the existence of magnetic fields (cf . IBorkowski et aLJ 
19901) . 

Since no self-consistent radiation-hydrodynamics computa- 
tions including heat conduction have been performed to datqj 
and urged by the fact that the new observations which became 
available by the Chandra and XMM-Newton satellites lack a 
convincing interpretation by detailed modeling, we decided to 
update our ID radiation-hydrodynamics code NEBEL by in- 
corporating a thermal cond uction module. We rec omputed then 
some of our sequences fromj Perinotto et aLJ (i2004 Paper I here- 
inafter) with heat conduction self-consistently included. Based 
on the new thermal structure of the bubble we modelled the X- 
ray and EUV emissions at selected positions along the evolution- 
ary sequences by using the well-documented CHIANTI code. 

We emphasize that in our study (i) wind and stellar evolution 
are consistently connected within the framework of the theory 
of radiation-driven winds, and (ii) no fit to the data by manipu- 
lating the models is done. It is our belief that only with such an 
approach it is possible to interprete the data appropriately and 
to deduce meaningful constraints regarding the relevant physics. 
Obviously, our approach is only applicable to objects which do 
not depart too much from sphericity, have closed shells where a 
bubble of shocked, hot wind gas is able to exist, and have central 
stars of normal, hydrogen-rich surface composition. 

We avoid to discuss in the present work objects with [WC] 
central stars for two reasons: (i) Their evolution is at present not 
known, and (ii) the physics of heat conduction in hydrogen-free 
and carbon-rich plasmas has still to be worked out. We note that 
the X-ray emission from a hydrogen-free plasma is expected to 
be stronger, making it more likely to detect extended X-ray emis- 
sion from PNe with [WC] central stars. The present number of 
positive detections is, however, too small as to make any defini- 
tive statements concerning the ratio of PNe X-ray sources with 
normal and [WC] central stars, respectively. 



The only exception is the central star of NGC 2392 whose wind 
speed is much too low for its position in the Hertzsprung-Russell dia- 
gram. 

^ JMellema & Frankl ( Il995t) computed the X-ray emission of their 2D 
hydrodynamics models, but heat conduction was not considered. 



Table 1. Elemental abundances, ei, used in the computations of 
our hydrodynamical models, in (logarithmic) number fractions 
relative to hydrogen, log e, = log(«i/«H) + 12. 



H 


He 


C 


N 


O 


Ne 


S 


CI 


Ar 


12.00 


11.04 


8.89 


8.39 


8.65 


8.01 


7.04 


5.32 


6.46 



The computational details of how the thermal conduction is 
treated are described in the Sect. |2l followed by a discussion of 
the resulting bubble structures in Sect. |3] Section |4] is devoted 
to the X-ray emission emerging from the conduction front and 
how it develops with time. In Sect. |5] we discuss extensively how 
our models compare with the existing observations. The paper 
conclud es with Sect. [6l Part of the results presented here can be 
found in ISchonberner et al.1 (l2006t) . 



2. The computations 

2. 1 . The hydrodynamical models 

The basic idea behind our modeling is to couple a spherical cir- 
cumstellar envelope, being the relic of the strong wind on the 
asymptotic giant branch (AGB), to a post-AGB star model of 
certain mass and to follow numerically the hydrodynamical evo- 
lution of the envelope across the Hertzsprung-Russell diagram 
towards the white-dwarf cooling path. This is achieved by apply- 
ing our ID radi ation-hydrodynamics code NEB EL as described 
in more detail in P erinotto et al.l (11998). The inner boundary con- 
dition in terms of density and velocity is provided by the variable 
central star wind. For the ionizing photon flux the star is assumed 
to radiate as a black body for each given eff'ective temperature. 

The radiat ion part of our hydrodyna mics code, CORONA, 
is described in lMarten & Sczcerbal (Il997h . We point out that this 
code is designed to compute ionization, recombination, radiative 
heating and line cooling fully time-dependently. At each volume 
element, the cooling function is composed of the contributions 
of all the ions considered and computed according to the actual 
plasma parameters. For each chemical element listed in Table[Tl 
up to 12 ionization stages are taken into account, amounting to a 
total o f 76 ion s. 

In ' Paper ll a very detailed description of how the nebular 
properties depend on the chosen initial envelope configurations 
and central-star models is given. For the present study we se- 
lected representative sequences and recomputed them with and 
without heat conduction included. The new simulations consider 
9 chemical elements instead of 6, with abundances as listed in 
Table [T] A compilation of the sequences investigated here is 
given in Tabl e El The 0.595 Mrp post-A GB model has been in- 
troduced in Schonberner et al.l ( 2005bl) . Sequence No. 10a is a 
variant of sequence No. 10 of iPaper a with the AGB mass-loss 
rate doubled. Note that the luminosities listed in Table |2] corre- 
spond to the early part of the post-AGB evolution. They decrease 
slowly during the evolution across the Hertzsprung-Russell dia- 
gram. 



Table 2. Hydrodynamical sequences of model planetary nebulae 
us ed in th is work. The s equence numbers correspon d to Table 1 
in'Paper I and Table 2 in ISchonberner et al.l (l2005bl) . but the ad- 
ditional notation 'HC and 'HC2' indicates that heat conduc- 
tion is included according to method 1 and 2, respectively (see 
Sect. 12.2b . The stellar luminosities refer to T^ff - 30 000 K, and 
the peak mass-loss rate of the AGB hydrodynamical simulation 
(6, 6a) is about 1 x lO"'* MQyr"^ Type indicates the structure 
adopted for the AGB envelope: 'A' means constant mass loss 
rate, p ^ r^^, 'B' means structure from hydrodynamical simula- 
tion (see Paper I for details). 



No. M L Magb 

(Mo) (Lq) (Moyr-') 



(kms" 



Type 



22 


0.565 


3981 


3 X 10-^ 


10 


A 


22-HC2 


0.565 


3981 


3 X 10-5 


10 


A 


6a 


0.595 


5 593 


Hydro, sim. 


= 12 


C 


6a-HC 


0.595 


5 593 


Hydro, sim. 


= 12 


C 


6a-HC2 


0.595 


5 593 


Hydro, sim. 


= 12 


c 


6 


0.605 


6280 


Hydro, sim. 


= 12 


c 


6-HC 


0.605 


6 280 


Hydro, sim. 


= 12 


c 


6-HC2 


0.605 


6 280 


Hydro, sim. 


= 12 


c 


10 


0.696 


11615 


1 X lO-* 


15 


A 


10-HC2 


0.696 


11615 


1 X lO-"* 


15 


A 


lOa-HC 


0.696 


11615 


2 X 10-^ 


15 


A 


10a-HC2 


0.696 


11615 


2 X lO--* 


15 


A 



FollowinglSEitzeiJ d 19621) and lCowie & McKed (Il977l) . flie elec- 
tron mean free path A is a function of electron temperature, T^, 
and electron number density, n^, and can be written as 

A = 2.625 X 10^ TllnJXnK [cm], (2) 

where the Coulomb Logarithm, In A, can be approximated as 



lnA = 



9.425 + 3/2 InTe - 1/2 Inwe, Te < 4.2 x 10^ K, 



22.37 ■ 



In Te - 1 /2 In n^, T^> 4.2 x 10= K, 



(3) 



for a pure hydrogen plasma. The diffusion coefficient D is then 
given by 



D = 7.04 X 10"'' A «e tI'^ [ergs-' K"' cm"']. 



(4) 



At high Te and low n^, the electron mean free path A becomes 
very large according to Eq. ^ (actually it can become much 
larger than the dimensions of the 'hot bubble'), and the difiiision 
approximation is no longer valid. Rather, the heat flux cannot 
exceed the saturation limit 



^sat= 1.72x10-" Tl'^, 



[erg cm s ] . 



(5) 



This is an approximate upper limit expressing the fact that the 
heat flux cannot be larger than the heat content, 3/2 n^kTs, times 
a characteristic electron transpo rt velocity, t^char- For a more de- 
tailed derivation of Eq. (HJ see lCowie & McKeel (Il977h . 



2.2. The treatment of electron heat conduction 

2.2.1. Physical description 

Electron heat conduction is described as a difiiision process, with 
the heat flux q given by 



q^-D VT,. 



(1) 



2.2.2. Numerical Treatment 

Each time step At of the hydrodynamical simulations is divided 
into 3 successive steps for updating the state vector Q (operator 
splitting): 

(a) Q(t) => advection ^ Qi(t + At); 

(b) Qi(t + At) ^ heat conduction => Qzit + At) 



(1st energy update at constant mass density p); 
(c) Qi(t + Af) => radiation ^ Q (t + At) 
(2nd energy update at constant mass density p). 

Ionization is frozen during step (b), being updated subsequently 
in step (c). For step (b), we solve the diffusion equation in spher- 
ical coordinates. 



dt 



■ PCo 



'df 



1 -5 , z 
-r —Ir"^ D 



dr 



dr 



(6) 



(E: internal energy per unit volume, c„: specific heat at constant 
volume, r. radial coordinate) with a fully implicit, standard nu- 
merical method. No additional constraints are imposed on the 
time step by this implicit energy update. 

The diffusion coefficient D is evaluated at the cell bound- 
aries from the local physical conditions according to Eq. (|4]i, 
with A replaced by A (see below) to take into account the effect 
of saturation. In order to take care of the fact that the diffusion 
coefficient itself depends on temperature, we have adopted the 
following procedure to obtain the temperature update from ini- 
tial temperature T[ to final temperature T2 due to electron heat 
conduction acting over the time interval Af : (i) we calculate the 
initial diffusion coefficient Do = D{T[y, (ii) starting from Ti, 
we solve the diffusion equation Eq. Q for time step Af/2 with 
D = Dq taken to be constant over time, resulting in intermedi- 
ate temperature ri/2; (iii) we calculate the intermediate diffusion 
coefficient D1/2 = D(Tii2); (iv) starting again from T\, we solve 
the diffusion equation for the full time step Af with D - D1/2 
taken to be constant over time, resulting in the final tempera- 
ture state T2. This procedure turned out to be perfectly adequate: 
comparison with a test calculation using 10 intermediate time 
steps instead of 1 gave practically identical results. 

In the framework of the diffusion approximation, Eq. ([T]i, sat- 
uration effects can be crudely accounted for by limiting the mean 
free path A. However, there is some arbitrariness involved in this 
approach, and it is not clear which recipe gives the most realistic 
results. We have tried two extreme methods: 

In method 1, we Umit the electron mean free path to a frac- 
tion of the spatial resolution of the numerical grid, i.e. we calcu- 
late A as 



Ai = min{/ x Ar, 2.625 x 10^ TljnJ In A), 



(7) 



where Ar is the local spacing of the radial grid. Choosing the 
constant / - 0.244 ensures that the conductive heat flux can 
never exceed the saturation heat flux given by Eq. (5): In regions 
where the limiter is active, i.e. where A\ - 0.244 ■ Ar, we obtain 
q - 9sat ■ Ar VT^e jT^ x q^^^ ■ ATf^/Tf., where AT^e is the tem- 
perature difference between two adjacent grid points. Hence, the 
limiting flux can only be reached at the sharp edge of a hot re- 
gion where lAT^el ~ T^. A drawback of this method is that Ai, 
and hence the results, depend explicitly on the numerical resolu- 
tion Ar. Sequences computed with this approach are denoted in 
TableE] with label 'HC. 

In method 2, we limit the electron mean free path according 
to the interpolation formula 



1 

I2 



1 



2.625x105 T^/nJ In A 



H- 4.105 



IVFe 



(8) 



As with the first method, the conductive heat flux cannot exceed 
the saturation heat flux given by Eq. (5): In regions where A <K 
Ar, the first term on the right hand side of Eq. ^ dominates, and 
A2 - A, i.e. the limiter is inactive. When A s> Ar, the second 



term on the right hand side of Eq. dSJ dominates, so A2 - 0.244 ■ 
Te/ IVFel, and q - q^^f. Sequences with this treatment of thermal 
conduction are indicated by the label 'HC2'(see Table|2]). 

The two methods represent the extreme cases. Method 2 es- 
sentially gives the saturation flux Eq. (|5]l wherever the unlimited 
flux would exceed this value, while method 1 generally yields 
much smaller fluxes, reaching the saturation flux only under ex- 
treme conditions. 

As a final remark we note that we have not attempted to cor- 
rect the calculation of the diffusion coefficient for the fact that 
the actual chemical composition of the nebular matter is not pure 
hydrogen. 

2.3. The CHIANTI code and the computation of the X-ray 
emission 

For obtaining the X-ray emission of our model PNe it is nec- 
essary to compute a synthetic optically thin X-ray spectrum for 
each radial shell of the model, since the different values for Tg, 
He and Mp at each shell result in unique spectral characteristics. 
The spectra were calc ulated using the CHIANTI software pack- 
age (iDere et al.lll997h . which has been used extensively by the 
astrophysical and solar communities. CHIANTI consists of an 
up-to-date set of atomic data for a large number of ions of as- 
trophysical interest, and also includes a number of ancillary data 
and a suite of useful routines. 

Our sy nthetic s pectra w ere co mputed with version 5.1 of 
CHIANTI (iLandi & Phillips l 2005l) . which includes the most re- 
cent atomic data. Only radial shells with Tg > lO'' K were con- 
sidered, since the contribution of cooler layers is negligible in 
the X-ray range. With the input of T^., n^ and n^ for a given 
volume element, individual spectra are then synthesized, includ- 
ing the contribution due to lines and various continua (free-free, 
free-bound an d two-photon co ntinuum). We have used the ion 
fractions from iMazzotta et al.l ( Il998.) . under the assumption of 
ionization equilibrium. The elemental abundances used are ei- 
ther those listed in Table[T| or solar values for the elements not 
considered in our hydrodynamics simulations. 

Each individual spectrum actually represents the emission 
(in erg s"' cm"^ A"') from a unit volume. In order to determine 
the total X-ray emission, the averaged spectrum, and the bright- 
ness distribution in the plane of the sky, we must perform ap- 
propriate integrations over our spectra, taking the (spherical) ge- 
ometry of the model into account. The multiplication of an in- 
dividual spectrum with the volume of the corresponding radial 
shell gives the total emission from that shell, and the summation 
over all shells emitting in X-rays then yields the total spectrum. 
An integration over a certain wavelength range then results in a 
well-defined total X-ray luminosity, Lx (in erg s ' ). 

For determining the brightness distribution in the plane of 
sky (in erg cm"^ s"' sr"'), the integration is performed along the 
line of sight (perpendicular to the plane of sky) for a series of 
impact parameters (up to about 400). 

The computation of the complete X-ray spectrum and lumi- 
nosity is rather time consuming and thus only performed for 
selected models along an evolutionary sequence. Such a post- 
facto computation of the X-ray emission is somewhat inconsis- 
tent with the hydrodynamics since at least part of the energy lost 
from the bubble by X-rays, viz. the energy loss due to ions not 
considered in our simulations (cf. Table[T] and Sect. 12.11 ). is not 
included in its energy budget. However, since the X-ray lumi- 
nosity computed by means of the CHIANTI code is only a very 
small fraction of the total radiation losses from the hot bubble 
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Fig. 1. The fractions of the last 4 ionization stages of oxygen 
within a typical PN bubble selected from sequence No. 6a-HC2. 
The model parameters are L = 5205 L0 and reff = 71667 K 
at t-5642 yr. Heat conduction is considered according to 
method 2. The central star is at the origin, the (reverse) 
wind shock at r-3.2 x 10'^ cm, and the conduction front at 
1.22 X 10'^ cm. Top: predictions from CHIANTI. Note that the 
CHIANTI code is only applicable to the region embraced by 
the wind shock and the conduction front. Bottom: predictions 
from our time-dependent computations with NEBEL. Beyond 
the conduction front, i.e. in the 'cool' nebular gas, the fractions 
of the displayed ions drop rapidly to virtually zero. 



as already provided by the NEBEL code (see Sect. l3.3l Fig.|7]i, 
a somewhat incomplete consideration of the X-ray loss in the 
bubble's energy budget has absolutely no consequences for the 
dynamics of the whole system. 

We checked also whether the ionization structures predicted 
by our hydrodynamics code and by CHIANTI are consistent 
with each other. Figure [T] gives an example for the ionization 
structure of oxygen within a typical bubble and compares the 
predictions of the CHIANTI code (top panel) with our time- 
dependent NEBEL/CORONA code (bottom panel). Note that the 
ionization fractions computed by CHIANTI are equilibrium val- 
ues and a function of the electron temperature only. In contrast, 
the ion densities computed by NEBEL/CORONA are the result 
of solving explicitly the time-dependent rate equations for the 
local temperatures and densities, and accounting for advection. 
Hence, the NEBEL results are expected to approach the equi- 
librium solution provided by the CHIANTI code only for suffi- 
ciently high electron density and low flow velocity. 

The bottom panel of Fig. [T] illustrates nicely how the ion- 
ization of oxygen changes while the gas passes through the 



wind shock: The main ionization stages of oxygen in the freely 
streaming wind are those of O^^ and O^^. After passing through 
the wind shock, the ionization switches to O^^, O^^ and O^^ be- 
cause of the very large post-shock temperatures. Also close to 
the conduction front where heated ma tter streams inwards (rel- 
ative to the front, see iBorkowski et a l. 1990, Fig. 4), we see a 
small influence of advection (compare the distribution of O''^ 
around r = 1.2 x 10'^ cm in both panels of Fig.[TJ. 

The dense regions above 10'' K close to the conduction front, 
i.e. for r k, 10'^ cm in Fig. [T] contribute most to the bubble's 
X-ray emission, as we will see later There we have reason- 
able agreement between the CHIANTI and NEBEL predictions. 
The remaining diff'erences are likely due to differences in the 
atomic data used. Considering the other uncertainties involved 
in this investigation, like, e. g., the wind model, the assumption 
of sphericity, and, last but not least, the observational data, we 
think that our approach of computing the X-ray emission post- 
facto from the nebular models by using the CHIANTI code is 
well justified. 



3. Heat conduction and bubble structure 

In this section we discuss in detail the results following from 
our numerical concept introduced in Sect. l2.2l and compare them 
with our previous simulations without heat conduction. For this 
purpose we used the sequences Nos. 6a, 6a-HC and 6a-HC2, all 
of which are based on the 0.595 M0 post-AGB model with the 
same initial envelope (cf. Table |2]l but with a different treatment 
of heat conduction. 



3.1. The wind model 

Figure |2] illustrates the evolutionary properties of the cen- 
tral star and its wind in terms of post-AGB time and stel- 
lar effective temperature. The relevant quantity for powering 
any X-ray emission is the mechanical luminosity of the stel- 
lar wind, Lwind =A^wind y^in(j/2. It is important to emphasize 
that, according to the theory of radiation-driven winds for stan- 
dard hydrogen-rich ch emical composition in the formulation of 
iPauldrach et al.l (Il988h . the mass-loss rate and the wind speed 
depend on the stellar parameters (mass, luminosity, effective 
temperature). Based on these wind prescriptions, the mechan- 
ical energy transported by the wind increases during the evo- 
lution across the Hertzsprung-Russell diagram, simply because 
the slowly decreasing mass-loss rate is over-compensated by 
the increasing wind speed (Fig. [21 upper right). However, when 
the hydrogen shell becomes exhausted, the mass loss rate drops 
sharply in line with the stellar bolometric luminosity, causing 
also the mechanical wind power to drop considerably since the 
wind speed remains now virtually constant at its maximum value 
of about 10 000 km s-^ 

In any case, the mechanical power remains always rather 
small and, in this particular case, does not exceed 1 % of the stel- 
lar photon luminosity (Fig.|2l lower left). According to the wind 
model used in this work, the maximum of the mechanical power 
occurs close to maximum stellar temperature. Only very little 
mass is carried away by the wind during the whole transition to 
the white-dwarf domain, viz. *3 x 10 "* Mq (lower right panel), 
which may be compared with the typical PN mass of a tenth of 
a solar mass. 

We emphasize that most of this mass is already lost with 
low speed during the first 1000 years of transition to the PN 
stage! Only during this phase we have wind speeds as low as 
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Fig. 2. Top: evolutionary path of the 
0.595 Mq model in the Hertzsprung-Russell 
diagram with the post-AGB ages, t, indi- 
cated along the track (left), and the stellar 
bolometric (solid) and wind (dashed) lumi- 
nosities vs. age (right). Middle: mass-loss 
rate (left) and (terminal) wind velocity 
(r i^ht). We followed the recommendations 
of lPauldrach et alJ ( Il988h for the central-star 
wind (Teff > 25 OO P K), while we assumed 
a Reimers wind (i ReimersI Il975h during 
the tra nsition to the PN domain (cf. also 
iPaper B) . Bottom: stellar wind luminosity 
(power) vs. stellar effective temperature 
(left), again with the post-AGB ages indi- 
cated, and the total mass lost by the wind 
during the post-AGB evolution (right). 
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Fig. 3. Stellar wind luminosities vs. effective temperatures for 
the four mass sequences considered in this work. The observed 
wind luminosities of PNe with diffuse X-ray emission, taken 
from the compilation in Table[3] are shown as 'star' symbols. 



a few 100 kms ', i.e. low enough to provide post-shock tem- 
peratures of the order of 10^ K. These mass-loss parameters are 
typical for thos e of the 'early wind' and are here based on the 
iReimersI (Il975h prescriptions (see caption of Fig.[2]l. The 'early- 
wind' phase is included consistently in our simulations with the 
appropriate treatment of radiative cooling. The stellar mass lost 
during the following PN stage is only about 8 x 10"^ M0, but 
this material has a very high kinetic energy because of its large 
speed exceeding 1000 kms ', leading to post-shock tempera- 
tures in excess of 10^ K. 



We close this section on our wind model with a discussion 
about its relevance for real objects. For this purpose we com- 
pared the observed wind luminosities of PNe with diffuse X-ray 
emission used in Sect.|5]and listed in Table [3] with the predic- 
tions of our post-AGB models (see Fig. |3]l. There is only one 
apparent discrepancy between theory and observation: The wind 
power of NGC 2392 is about a factor of 5 below our predictions, 
a consequence of the exceptionally low wind speed measured 
for this particular object. The remaining objects from Table [3] 
(NGC 3242, NGC 6543, NGC 7009) have wind luminosities 
which are, on the average, only a factor 2 below the theoretical 
predictions. 

Given the large uncertainties of the mass-loss rate determina- 
tions, we do not consider a factor of two difference between ob- 
served and computed wind powers to be alarming. Moreover, the 
distances of NGC 3242, NGC 6543, and NGC 7009 used here 
provide stellar luminosities that are somewhat below those for a 
typical central star of 0.6 Mq which is about 5000 Lq (cf. Table 
|3]l. A coiTesponding increase of the distances would bring theory 
and observation in much closer agreement (L^ind '^ distance''^). 



3.2. The influence of thermal conduction 

The influence of thermal conduction is illustrated in Fig.|4]where 
the density and temperature structures of three models at about 
the same position along the evolutionary path shown in Fig.|2] 
are compared. As expected, thermal conduction across the bub- 
ble/PN interface has a profound impact on the density and ther- 
mal structure of the bubble. As a reference, the top panel of 
Fig. |4] shows the typical temperature / density structure of the 
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Fig. 4. Radial profiles of electron density (solid, left ordi- 
nate) and temperature (dashed, right ordinate) of three models 
taken from sequences No. 6a (top), No. 6a-HC (middle) and 
No. 6a-HC2 (bottom), respectively, at about the same positions 
along the stellar path shown in Fig.|2l The approximate stel- 
lar parameters are L ^ 5200 Lq, Jeff - 71 600 K at f :^ 5640 yr. 
The central star is at the origin, and the (reverse) wind shock at 
r ^ 3 X 10'^ cm. The shocked wind gas, i.e. the bubble, is (in 
all three cases) between the wind shock and the contact sur- 
face/conduction front (r ^ 1.3 x 10'^ cm). The PN proper is 
bounded by the contact surface/conduction front and an outer 
shock at r ^ 4.4 x 10'^ cm. The PN is surrounded by the ion- 
ized AGB wind whose radial density profile reflects the mass- 
loss history of the late AGB evolution (ISchonberner et al]ll997t 
ISteffenetal.lll998h . 



shocked wind if thermal conduction is ignored. The gas is very 



hot (Te ^ 8 X 10^ K) and tenuous (Wg ~ 0.5 cm""*), and these 
properties do not change much with radius. Such a bubble struc- 
ture is in sharp contrast to the cases where heat conduction is 
explicitly considered, as is evident from the middle and bottom 
panels of Fig.|4] 

Akeady the treatment of heat conduction according to 
method 1 leads to a completely different bubble structure: a sig- 
nificant temperature gradient is established, whereby the region 
close to the conduction front reaches temperatures as low as 
10^ K (Fig.m middle panel). The bubble remains virtually iso- 
baric (because of the high sound speed), and an increased mat- 
ter density must compensate the temperature decrease. This ad- 
ditional mass is provided by a (subsonic) flow o f nebular gas 
heated ('evaporated') at the conduction front (cf. IWeaver et al.l 
Il977h . Note that the physical conditions immediately behind the 
(reversed) wind shock are virtually not influenced by heat con- 
duction, i.e. the shock remains fully adiabatic. 

Method 2 provides more efficient heat conduction, and hence 
the thermal structure of the bubble is more homogeneous with 
smaller temperature and density contrasts between wind shock 
and conduction front (Fig. |4] bottom panel). The bubble gas is 
not hotter than about 5 x 10^ K (except in a thin region im- 
mediately behind the wind shock), and the minimum (electron) 
density is about 10 cm"-'. Because thermal conduction carries 
energy away from the wind shock, the latter is now not adia- 
batic anymore: the density jump across the shock is nearly two 
orders-of-magnitude, in sharp contrast to the heat conduction 
treatment based on method 1 and the case without heat conduc- 
tion in which the densities increase only by a factor of four (see 
Fig.m top and middle panel). We note also that the bubble struc- 
tures seen in Fig. |3] are typical ones which do not change much 
along the main evolution across the HR diagram. The concept 
of an ideal contact discontinuity with no mass transfer from one 
side to the other does not hold anymore: The outer edge of the 
bubble (= inner edge of the nebula) is now defined by the heat 
conduction front. 

Figure|4]indicates also that the dynamics of the whole system 
remains practically unaffected by the physical treatment of the 
bubble gas: in all three cases considered, the bubble sizes and 
the nebular structures are virtually identical (see also Fig.fTOll. 

Heat conduction does not change the total energy budget 
of the bubble but only transfers thermal energy across the bub- 
ble from the wind shock towards the contact surface/conduction 
front, where it is used to heat and 'evaporate' nebular gas. The 
latter flows inwards relative to the co nduction front and re mains 
thereby inside the bubble (see e.g. I Weaver et al.l [1977). Line 
cooling can change the bubble's energy content substantially 
(see Sect. l3.3l l. However, radiative cooling is confined to a thin 
surface layer of the hot bubble (10^ ... 10^ K) whose tempera- 
ture and density structure is not much affected by thermal con- 
duction. Radiative losses are therefore only slightly enhanced in 
models with heat conduction (cf. Fig.|7]i. Although the amount of 
mass 'evaporated' from the cold nebular gas exceeds by far the 
mass injected into the bubble by the stellar wind, it can be totally 
neglected for the mass budget of the PN proper (see Sect. 13. 3b . 

Figure|5]illustrates the typical bubble structures for two mod- 
els with different central stars and hence different wind powers 
and time scales of evolution. The model with the massive cen- 
tral star (0.696 M0) has a very compact and dense bubble due 
to its fast evolution and powerful stellar wind. The other ex- 
treme occurs for the model with the least massive central star, 
(0.565 M0), which evolves most slowly and has only a rather 
modest stellar wind: the bubble is extended and relatively ten- 
uous. Also the mean bubble temperatures reflect the different 
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Fig. 5. Radial profiles of electron density (solid, left ordinate) 
and temperature (dashed, right ordinate) of two models with dif- 
ferent central-star masses. Top: 0.696 M© (sequence No. 10- 
HC2) at f = 301 yr and with L = 1 1 392 Lq, T^s = 70 405 K. The 
large drop of electron density and temperature at r ^ 2.5 x 
lO'^ cm indicates the position of the ionization front. The re- 
verse wind shock is at r = 4.5 x lO'^ cm, the conduction front at 
r= 1 X 10'^ cm. Bottom: 0.565 M© (sequence No. 22-HC2) at 
f = 12 107 yr with L = 3455 L©, T^s = 69 906 K. The wind shock 
is at 1 X 10'^ cm, the conduction front at r = 4.5 x 10'^ cm. Note 
the different ranges the electron density. 



wind properties: the 0.696 Mq model has the largest bubble tem- 
peratures of our sequences because the mechanical energy input 
by the the stellar wind is largest and the bubble size is smallest, 
as explained in more detail in Sect. l3.3l The case of 0.595 M© 
shown in the bottom panel of Fig. His intermediate. 

3.3. The evolution of the bubble masses 

The evolution of the bubble's mass with time is shown in 
Fig. |6] for the two heat conduction cases considered here. One 
sees that the 'evaporated' nebular matter soon exceeds the mat- 
ter blown into the bubble by the stellar wind. After 10000 
years of evolution across the Hertzsprung-Russell diagram, the 
mass contained in the bubble reaches a maximum which de- 
pends on the treatment of thermal conduction: ^8 x 10 "* Mq 
(HC) and ^ 15 x 10 "^ Mg (HC2), respectively, which is still 
negligible compared to a typical PN shell mass of about 
0.1 M0. Obviously, the reduced thermal conduction efficiency 
of method 1 results in a lower evaporation rate of cold gas rela- 
tive to that of method 2. In either case, the whole bubble mass is 
virtually confined in a narrow outer shell. 
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Fig. 6. Evolution of the bubble mass with time for the sequences 
Nos. 6a-HC and 6a-HC2. Gas hotter than 10^ K is considered to 
belong to the bubble. The scatter is numerical noise. The dashed 
line indicates the mass blown into the bubble by the central-star 
wind, counted from t - 1200 yr. The dotted line gives a theoret- 
ical upper limit of the 'evaporated' mass based on Eq. ( fT6b . 

Note that most of the mass blown-ofFby the central-star wind 
during the post-AGB evolution is not contained in the hot bubble 
but in the inner part of the nebula. The (reverse) wind shock 
forms only when the wind speed exceeds a certain value, which 
happens at about 1200 yr after departure from the tip of the AGB 
in the case of a 0.595 Mq central-star model (for the zero point 
see Fig.|2]i. Only about 1 x 10""^ Mq are blown into the bubble 
during the remaining part of evolution (Fig.|6]l. 

The general time dependence of the bubble mass as shown 
in Fig. [6l can be in t erp reted in the following way. According 
to B orkowski et al.l (Il990.) . one can distinguish three different 
phases for thermal conduction fronts: (i) the 'evaporation' phase 
in which the front advances relative to the initial interface and 
heats cool gas; (ii) the quasi-static phase in which the front 
stalls because heating and 'evaporation' balance radiative cool- 
ing and 'condensation'; (iii) the 'condensation' phase in which 
the front recedes relative to the previous position because cool- 
ing/'condensation' dominates. 

As long as the wind power and hence the energy input into 
the bubble increases steadily, the conduction front is obviously 
in the 'evaporation' phase, and cooler nebular gas is heated and 
added to the bubble. This phase ends when the stellar lumi- 
nosity, and hence also the wind power, drops because hydro- 
gen burning ceases. This occurs in our simulations after about 
10000 years of post-AGB evolution for the 0.595 Mq model 
shown in Fig.|6](see also bottom left panel of Fig.|2]). For the fol- 
lowing a; 3500 years radiation cooling dominates, and the front 
becomes 'condensation' dominated, leading to a small reduction 
of the bubble mass (Fig. |6]l and the rate of expansion. After the 
rapid fading of the central star has stopped at about 12 000 years, 
the wind is blowing with nearly constant albeit lower strength 
(Fig.|2]i. Due to the continued expansion, densities decrease and 
radiative losses becomes less important; the conduction front 
turns slowly back into the 'evaporating' stage, and the bubble 
mass starts to increase again, but at a r educe d rate (Fig.|6]l. 

For spherical cases IWeaver et al.l (1 19771) estimated analyti- 
cally an evaporation rate of 



Mhb = Ci<rhb)'/'</(-Rhb- 
- Ci (J\iif^^R^h since 



Rs «J ^hb , 



(9) 



where Ci a; 4. 13x10 ''^ and(rhb), Rhb, and/?s are the mean tem- 
perature of the hot bubble, its radial size, and the position of the 
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Fig. 7. The ratio Qrad/^wind as a function of time for models 6a 
(thick) and 6a-HC2 (thin). Radiative losses are confined to a thin 
layer at the outer surface of the hot bubble; they are slightly 
larger in the presence of thermal conduction. 



wind shock, respectively. The central star defines the origin, and 
the wind shock at R^^ constitutes the inner boundary of the bub- 
ble. It must be noted, however, that Eq. (|9]l assumes that radiative 
cooling of the bubble is negligible, a nd thus provides only an up- 
per limit of the evaporation rate (cf. IWeaver et al.iri977l Eq. 61 
therein). 

Equation (|9]l allows us to compute the (equilibrium) evap- 
oration rate of a bubble with given temperature and radius. 
However, what we really want to know here is how the mean 
temperature of the hot bubble, and hence its evaporation rate, 
depend on the power of the stellar wind, Lwind - A^windf^injj/2. 
To answer this question, we consider the energy balance of the 
hot bubble, first without heat conduction. In this case, the ki- 
netic energy of the stellar wind is used to increase the thermal 
and kinetic engery of the bubble, to cover the work done by the 
expansion of the bubble, and to compensate possible radiative 
losses from the bubble: 



-^wind 



Eth + ^kin + Wexp + 2i 



irad, 



(10) 



where giad represents the net radiative cooling rate integrated 
over the volume of the bubble. Denoting by /th the fraction of 
the wind power that is converted into thermal energy, we can 
write 



3 k 

/th iwind - ^th - A^wind T " 7^0 , 
2 fl 



(11) 



where Mwind is the stellar mass loss rate and Tq the tempera- 
ture of the shock-heated bubble. For a strong adiabatic shock, 
/th=9/16 and Tq is given by the relation 



3 yU 2 



To ~ 1.4x10^ 



^wind 



1000 km/s 



[K] .(12) 



Here Uwind is the velocity of the stellar wind, and fi is the average 
mass of the gas particles. The numerical factor in Eq. (fT2l) was 
obtained assuming jj = 0.6 niu for fully ionized plasma of solar 
composition. According to Eq. (fT2l i. the temperature of the hot 
bubble depends only on the wind velocity. We note that the bub- 
ble temperatures of the two models shown in Fig. |5] cannot be 
understood with this relation: the stellar wind of the more mas- 
sive central star (top) is about two times slower than that of the 
less massive one (bottom), but its bubble temperature is about 3 
times higher! Obviously, thermal conduction makes a difference. 



In the presence of heat conduction, the mass of the of bubble 
increases both due to the mass loss through the stellar wind with 
rate Mwind, and due to 'evaporation' through thermal conduction 
with rate Mhb- Under these circumstances, the wind must supply 
the power to heat and expand both the stellar wind entering the 
bubble through the inner shock and the mass being added at the 
surface of the bubble by 'evaporation'. The mean temperature of 
the bubble, (Thb), is then given by the relation 



fthL. 



'wind 



3k_ 
3 k 



(M„i„d + Mhb) <rhb> 



= - - (Mwind (Tub) + Ci (rhb)'^' Rkb) , (13) 

where we have used Eq. ^. Equation (T3[ can be used to es- 
timate the mean temperature of the hot bubble for given wind 
power, wind mass loss rate, and bubble size. In turn, the evapo- 
ration rate follows from Eq. ^. 

If Mwind ■*: Mhb, which is a valid approximation during the 
main part of the evolution (once the bubble is well established 
until the central star begins to fade, see Figs.|6]and[8]), we obtain 



<rhb)^^^ = 



fthLv. 



(I^Ci^hb) 



or (rhb)^24 



'^wind 
^hb 



2/7 



(14) 



According to Eq. (fT4l i. the temperature of the hot bubble de- 
pends on the wind power and the bubble radius. This relation 
explains the behavior of the two models shown in Fig.|5j the 3 
times higher wind power of the more massive central star, and 
its 10 times smaller bubble size imply a mean bubble tempera- 
ture which is about 2.6 times higher that that of the less massive 
model. 

Combining Eqs. (|9| and (fT4l i. we can write the evaporation 
rate as 
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The dotted curves shown in Figs.|6] and |8] were obtained by in- 
tegrating Eq. ( fTSI l over time, given Lwind(0 and Rhbit) for the re- 
spective model sequence. Note again that Eq. ( fT6t gives an upper 
limit of the evaporation rate, since both /th and C] will be smaller 
than assumed here if radiation losses are important. 

In our model simulations, 'evaporation' due to thermal con- 
duction is significantly less efficient than predicted by Eq. ( fTSI l. 
As shown in Fig.|6] the mass of the hot bubble measured in the 
numerical model is a factor 3 to 5 lower than suggested by the 
analytical estimate. In-depth investigations have revealed that 
this discrepancy is caused by substantial radiative energy losses 
at the conduction front / contact discontinuity. As illustrated in 
Fig-El these losses are most severe in the early phases, when a 
large fraction of the wind power is radiated away. This explains 
also why the bubble mass grows very slowly at the beginning: 
a considerable part of the mass added to the bubble by the stel- 
lar wind cools efficiently and 'condensates' outside the bubble. 
Later, 'evaporation' becomes more efficient and the mass of the 
bubble grows much faster than implied by the wind moss loss 
rate. But even at f = 8000 yr, 50% of the wind power is still 
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Fig. 8. The bubble masses vs. time for two sequences with a 
0.696 Mq central star (Nos. 10-HC2 and 10a-HC2) (top) and 
the 0.565 Mq sequence (No. 22-HC2) (bottom). Again, only 
gas hotter than 10^ K is considered to belong to the bubbles. 
The dashed line indicates the masses blown into the bubble 
by the central-star wind, counted from f = 75 yr (0.696 Mq) 
and t - 2600 yr (0.565 M0), respectively. The dotted line gives 
the theoretical upper limit of the 'evaporated' mass based on 
Eq.©. 



lost as radiation. It is interesting to note that the radiative cool- 
ing at the surface of the hot bubble is hardly altered by thermal 
conduction (see Fig.|7ll. 

For completeness we show in Fig. |8]the evolution of bub- 
ble masses for two sequences with a more massive (0.696 Mq) 
and for one sequence with a less massive central star. The two 
sequences Nos. 10-HC2 and 10a-HC2 with the massive central- 
star model (top panel) differ only by their AGB mass-loss rates 
of 1 X 10""* and 2 x lO""* Mq yr"', respectively (see Table|2]l. 

For the 0.696 M0 cases, the time evolution is similar as 
already discussed above for 0.595 Mq. A phase of rapidly 
increasing bubble mass is followed by a substantial decrease 
when the heat supply by the fading central-star wind drops. 
The 'condensation' of bubble mass is more severe in the 
0.696 Mq cases because the bubbles are small and rather 
dense (wg ~ 200 . . . 500 cm"-'), affording a more efficient radia- 
tive cooling. The discrepancy between the evaporation rate ob- 
tained from the simulation and the prediction of Eq. (fT6] l. respec- 
tively, is even more severe than in the case of the central star with 
0.595 Mq. 

Since the massive central-star models evolve on very short 
time scales, their bubbles are the smallest in size and mass of 
all sequences computed. The slightly smaller bubble masses 
found for sequence No. 10a-HC2 (as compared to those of se- 



quence No. 10-HC2) are due to their somewhat smaller sizes, 
a consequence of the denser bubble environment caused by the 
larger AGB mass-loss rate, M^gb = 2 x 10"'* Mq yr"' instead of 
Magb = 1x10-4 Mgyr"'. 

In the case of the low-mass central star (0.565 Mq), the wind 
power drops slowly after 20 000 years of post-AGB evolution. 
As a result, the mass of the hot bubble ceases to increase. But 
in contrast to the simulations with more massive central stars, 
there is no real 'condensation' phase. Rather, the mass of the hot 
bubble remains constant with time, obviously because radiative 
cooling is inefficient due to the low plasma densities. Hence, we 
find the best agreement between analytical and numerical evap- 
oration rate for this sequence. 

We finally note that 'evaporation' due to heat conduction 
(HC2) increases the bubble mass by as much as a factor of 10, 
relative to the wind-blown mass, for the 0.696 Mq cases, and 
by a factor of about 40 for the model with M=0.565 Mq near 
the end of the evolution (t > 20 000 yr), when 'evaporation' and 
'condensation' balance and the bubble mass is essentially con- 
stant (see Fig.|8]l. The case of 0.595 Mq is intermediate. 



4. The X-ray emission 

4.1. Spectra 

We computed detailed X-ray spectra, as described in Sect. 12.31 
for the wavelength interval 5-100 A, corresponding to 2.5- 
0.125 keV. Examples are shown in Fig.|9]for the three models 
from Fig. HI The X-ray emission consists of a continuum and 
numerous strong lines of highly ionized species. The models 
with thermal conduction included (middle and bottom panel of 
Fig.|9]l have much larger fluxes in the 5-100 A wavelength band 
than the model without conduction. The total X-ray luminosi- 
ties (in the same spectral band) are 2.4 x 10-" erg s"' (HC) and 

4.2 X 10^' ergs"' (HC2), coiTesponding to 1.2 x 10"^Lstai- (or 

6.3 X 10"^ Lq) and 2.1 x 10"''Lstai- (or 1.1 x 10"^ Lq), respec- 
tively. Without thermal conduction, the X-ray luminosity is only 
2.7 X 10^" erg s"' in the same wavelength band. In any case, these 
luminosities are only minute fractions of the mechanical energy 
carried away by the central-star wind, which is about 3.4 Lq for 
the models used in Fig.|9](cf. Fig.|2]). 

We note that our models with thermal conduction included 
provide X-ray luminosities which ar e smaller than those com- 
puted by IZhekov & Perinottol (Il996h (by roughly a factor 10). 
However, radiative line cooling is not included in their analyti- 
cal models, leading to a severe overestimation of the evaporation 
rate (see Sect. l3.3l ). Moreover, the simplified assumptions that 
only wind interaction is responsible for the compression of the 
AGB wind, ignoring the effects of photoionization, result in an 
unreahstic expansion behavior of the bubble. 

We conclude th a t an a pproach like the one conducted by 
IZhekov & Perinottol (Il996h is likely too simple for a reliable 
computation of the X-ray emission from planetary nebulae. 

We emphasize in this context that the X-ray flux depends 
also on the chemical composition of the bubble gas. A de- 
tailed investigation of how the X-ray flux depends on the chemi- 
cal composition, especially for cases with hydrogen-poor stellar 
winds, will be the subject of future work. 



4.2. Surface brightness distributions 

The X-ray intensity, integrated over the wavelength band of 5- 
100 A, is limb-brightened and reflects the fact that the regions 
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Fig. 9. The fluxes emitted in the wavelength band of 5-100 A 
(2.5-0.125 keV) computed from the bubbles of the models 
shown in Fig. |4] without heat conduction (top), and with heat 
conduction included: HC (middle) and HC2 (bottom). The fluxes 
given in the ordinate refer to an assumed distance of 1 kpc, and 
the contribution of regions with Tg < 10^ K is ignored. The in- 
sets give the luminosities in three diff'erent wavelength bands. 



with physical properties favorable for emitting X-rays are mainly 
the outer layers of the hot bubble (Fig. [TOl i. The absolute value 
of the intensity and its radial profile depend therefore sensitively 
on the way thermal conduction is treated. 

In the case without any thermal conduction, the X-rays 
are emitted only from the very thin interface between the hot 
shocked stellar wind and the PN gas, resulting in a ring-like, 
comparatively weak emission. (Fig.[TOl top panel). In the case of 
low-efficiency heat conduction (method 1) the X-ray emission is 
still rather strongly peaked towards the bubble's surface, with a 
center-to-limb variation of a factor 4 (middle panel), method 2 
results in a more homogeneously distributed X-ray emission 
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Fig. 10. Radial intensity profiles, integrated over the wavelength 
band of 5-100 A (2.5-0.125 keV), of the models from Fig. |9] 
without (top) and with heat conduction considered: HC (mid- 
dle) and HC2 (bottom). At impact parameter p-3.3 x 10'^ cm 
the inner edge of the bubble, i.e. the position of the (reverse) 
wind shock, is visible in the bottom panel. 



with a modest center-to-limb variation of only a factor 1.5. A 
larger part of the bubble is now contributing to the X-ray emis- 
sion (bottom panel). 

One notices that, although all the three models shown in 
Fig. [To] originated from the same initial configuration and are 
virtually of equal age, the sizes of their bubbles differ by a small 
amount: with thermal conduction included the bubble is a little 
bit smaller, rhb = 1.2 x 10'^ cm instead of rhb = 1.3 x 10'^ cm 
in the case without conduction^. This appears to be related to the 
fact that radiative energy losses from the bubble are somewhat 
larger for models including thermal conduction (see Fig.|7]l. This 



In the cases with heat conduction, the size of the bubble, ri,b, 
defined by the outer edge of the conduction front. 



additional cooling causes a slightly slower expansion of the bub- 
ble, but is too small to significantly change the dynamics of the 
whole nebula: The PN evolution is virtually not influenced by 
heat conduction (cf. Fig.|4]i. 

4.3. Time evolution 

An important criterion for the usefulness of thermal conduc- 
tion in explaining the X-ray emission from planetary nebulae 
is not only the predicted X-ray luminosity itself, but also how it 
evolves with time, i.e. as a function of the stellar parameters. In 
general, the X-ray luminosity is determined, for a given temper- 
ature, by the total volume emission measure of the X-ray emit- 
ting region, rjx - pxMx, with px being a characteristic mean 
density of the emitting volume and Mx its mass. For simplic- 
ity we set Px K MxIR\^. If Mx oc f " and /?hb °^ f *, we have 
Lx oc pxMx oc f"-^'' . 

Considering for the moment no thermal conduction, and as- 
suming a central-star wind with constant mass-loss rate and ve- 
locity, we have a - \ and b - \, where the latter expression 
holds if the bubble ex pands into an environment with p oc r^^ 
(lKoo&McKe3ll992l Eg. 3.1 therein). In this case, Lx ^ r' 
(see also IVolk & Kwoklll985L Eq. 33). This situation changes 
somewhat if the X-rays are only observed within a limited en- 
ergy range and if the wind is evolving in time, i.e. if the wind 
speed increases (cf. Fig. |2]i. Since the post-shock temperature 
scales with wind speed squared, the X-ray emission increases 
first, reaches a maxi mum at an appropria te wind velocity, and 
drops then rapidly (see Volk & Kwokl 19851 Figs. 1 1-13 therein). 
Note that in the case of increasing stellar wind power the bub- 
ble expansion is accelerated {b > 1), leading to a faster X-ray 
luminosity drop than in the case of a constant wind. 

We have already seen in Sect. [3] that thermal conduction 
across the bubble/PN interface causes the bubble mass to in- 
crease rapidly with time thanks to 'evaporation' from the dense 
nebular matter. For the case shown in Fig. |6] a - 4.5 and b ^2, 
leading to Lx oc p. Any variation of Tx, the characteristic tem- 
perature of the X-ray emitting volume, during the course of evo- 
lution is thereby neglected. 

This variation of the emission measure is the consequence 
of the adopted wind model in which wind power and speed in- 
crease during the evolution as shown in Fig. |2] Obviously the 
bubble mass increases faster by heat conduction than the density 
decreases by expansion. In contrast, the X-ray luminosity would 
increase only linearly with time if we assume Mhb '>^ Rhb (Eq.|9j 
and a linear expansion law, /?hb '^ t, instead. 

According to the approximate expression for the emission 
measure given above, it is also possible that Lx decreases with 
time, provided 2a - 3b becomes negative, or a < 1.5 b, at some 
point during the evolution. This happens for instance in our sim- 
ulations when the wind power decreases at the end of the evolu- 
tion (a becomes even negative, cf. Figs.|6]and[8]l. Another possi- 
bility is a weak wind generating only a small evaporation rate. 

The temporal evolution of the X-ray emission for a selection 
of our sequences listed in Table[2] is presented in Fig.[TT] We 
consider only models computed with heat conduction according 
to method 2 (HC2). All sequences are plotted from the early 
transition phase until the white-dwarf domain is reached. Note 
that the evolutionary time scales decrease rapidly with stellar 
mass; it takes only ^ 1000 yr for the 0.696 M© central star to 
exhaust its hydrogen-burning shell. 

At first glance it appears to be astonishing that, according 
to our simulations. X-ray luminosities increase systematically 
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Fig. 11. Development of the X-ray luminosity in the wavelength 
range 5 to 50 A (2.5-0.25 keV) as a function of time for the hy- 
drodynamical sequences indicated in the inset. Five sequences 
have heat conduction included (22-HC2, 6a-HC2, 6-HC2, 10- 
HC2, and 10a-HC2 of Table |2|, and three are without heat con- 
duction for comparison (6a, 6, and 10, dashed). The strong fluc- 
tuation of X-ray luminosities of these models reflect the large 
numerical uncertainties related to the thinness of the X-ray emit- 
ting region. The symbols indicate the times at which the X-ray 
emission was computed by means of the CHIANTI code (cf. 
Sect. 12.31 ). The two panels show Lx vs. time (top) and Lx/^wind 
vs. time (bottom), respectively. 



with central-star mass, despite the fact that the bubble masses 
decrease (cf. Sect 13.3b . This is, however, a consequence of the 
dependence of post-AGB evolution on central-star mass: Both 
the stellar wind power and evolutionary speed increase substan- 
tially with remnant mass. As a result, the hot bubbles around 
more massive objects remain smaller and more dense and have 
larger emission measures despite of their smaller masses as com- 
pared to those around less massive and more slowly evolving 
central stars. 

The trend of the X-ray emission with age is expected from 
the previous discussions. At early times, at the beginning of 
the sequences shown in Fig.[TT] the X-ray luminosities increase 
sharply with time, reflecting the increase of wind speed and the 
resulting increase of post-shock temperatures to values suitable 
for soft X-ray emission (top panel). Note that both types of mod- 
els, with and without heat conduction, give virtually the same 
X-ray emission during this phase, because the fraction of 'evap- 
orated' matter is still insignificant. This is typical for the 'early- 
wind' phase of evolution when heat conduction is still unimpor- 
tant. 



For the energy range considered here (0.25-2.5 keV), 
there appears to be an optimum velocity, ywind -400 kms"' 
(T^post-shock - 2.3 X 10^ K, see Eq.[T2]i for which the X-ray emis- 
sion of models without thermal conduction reaches a maximum. 
This happens at f ^ 95 yr for the 0.696 Mq sequence and 
t ^ 1800 yr for 0.595 Mq, respectively, when the central stars 
are still very cool {T^s < 25 000 K) and the planetary nebula for- 
mation is just beginning. The luminosities achieved (^ 10^"-10-" 
erg s"') are in the observed range, but the stellar parameters like 
temperature and wind speed are inconsistent with the observed 
values! Beyond this maximum, the X-ray emission decreases 
with time due to expansion and because the bubble becomes too 
hot in the absence of heat conduction. The memory of the 'early- 
wind' phase is lost. 

The situation changes completely if heat conduction is in- 
cluded in the simulations. The X-ray luminosity increases with 
time to levels up to 2 dex above the values achieved without con- 
duction. The trend with age is in line with the discussion above, 
i.e. an increase with time (or buble size) until the central star 
is rapidly fading (top panel of Fig.fTTTi. This occurs after about 
9000 years of evolution for 0.595 Mq, after about 6000 years for 
0.605 Mq, and after only 800 years for 0.696 Mq. 

The final evolution, once the fading of the PN nucleus has 
stopped, is interesting. We see from Fig. |6] that the heat con- 
duction turns back into the 'evaporation' stage (at t ^ 15 000 
yr for the M - 0.595 Mq example), albeit with much reduced 
efficiency. The 'evaporated' nebular mass is, however, not suf- 
ficient to compensate for the emission measure decrease due to 
expansion, and thus the X-ray luminosity decreases with time as 
well. 

Altogether we see that the numerical simulations confirm the 
qualitative estimates made above: During the horizontal part of 
the evolution across the HR diagram the high evaporation rate 
triggered by a powerful and accelerating stellar wind dominates 
and leads to an ever increasing X-ray power At the end of the 
evolution when the wind dies, the continued expansion reduces 
the emission measure for the X-rays. 

Figure[TT]demonstrates also the dependence of the X-ray lu- 
minosity on details of the final AGB mass-loss rate. First of all, 
the bubble of sequence No. 10-HC2 (Magb = 1 x lO^* Mq yr"') 
expands a Uttle bit faster than that of sequence No. 10a-HC2 
(Magb = 2 X 10""* MQyr"') because the latter has a denser en- 
velope. Consequently, they differ also somewhat in their X-ray 
luminosity. At the beginning of the evolution, the bubble of se- 
quence No. 10a-HC2 with its denser shell develops a slightly 
lower X-ray luminosity because its emission measure pxMx is 
lower: the smaller Mx is not adequately compensated for by a 
larger px- Later on, the radius difference between both sequences 
becomes larger, and the X-ray emission from sequence No. 10a- 
HC2 exceeds that of sequence No. 10-HC2 because of the large 
radius dependence of rjx- 

Our detailed hydro dynamical treatment does not confirm the 
scaling law derived bv IZhekov & Perinotta (1996) according to 
which Lx cc M"^^ for constant AGB wind speed and given age. 
This law predicts an X-ray luminosity for sequence No. 10a- 
HC2 (Magb = 2x10""* Mq yr"') which is larger than the X-ray 
luminosity of sequence No. 10-HC2 (Magb = 1 x 10"'*MQyr"') 
by 0.23 dex for the whole evolution (cf. also the previous discus- 
sion of this matter in Sect. l4.1b . Our dependence is smaller and 
more complicated. 

The bottom panel of Fig. fTTI illustrates in detail how the X- 
ray luminosity depends on the central-star's wind power. In gen- 
eral, a fraction of between 0.01 and 0.001 of the wind power is 



radiated away by X-rays between 5 and 50 A during the high- 
luminosity /high wind-power phase of evolution. These fractions 
become, of course, smaller (larger) if the wavelength band is re- 
duced (increased). The variations of Lx/iwind along a sequence 
reflect the response of thermal conduction to the stellar wind 
evolution: After the early maximum of X-ray emission the ratio 
Lx/^wind decreases first because it needs some time for the con- 
duction to 'evaporate' sufficient nebular matter to compensate 
for the expansion based reduction of the emission measure (^ 
pxMx)- Then the X-ray luminosity turns out to be fairly propor- 
tional to the wind power during the following evolution across 
the HR diagram. The sudden increase of Lx/^wind immediately 
before the stellar (wind) luminosity declines is due to the fact 
that 'condensation' occurs on a longer time scale than the stellar 
(wind) luminosity drop. 

Without heat conduction included, LxlL^ind becomes very 
small during the course of evolution, viz. ^ 10 ''-lO"^. This 
fact demonstrates that thermal conduction by electrons across 
the contact surface is an efficient means to convert wind power 
into X-ray radiation for the whole lifetime of a PN! 

There is an additional factor that influences the X-ray lumi- 
nosity: the spectral band used in the computations (and observa- 
tions as well). Here it is the low energy limit that matters, since 
energies higher than E^ 0.65 keV {A < 20 A) contribute very lit- 
tle to X-ray luminosity. Defining the X-ray band by 0.25 keV < 
E < 2.5 keV (5 A < /I < 50 A), more massive and luminous cen- 
tral stars show an apparently larger efficiency of converting wind 
power into X-ray power, as is clearly seen in Fig. [TT] Defining 
instead the X-ray band by 0.05 keV < £ < 2.5 keV (5 k < A < 
250 A), the ratio Lx/Lwind becomes roughly 0.01 for the whole 
range of central star masses. 



5. Comparison with observations 

In this section we will compare the results of our thermal con- 
duction simulations with existing observations from the X-ray 
satellites XMM-Newton and Chandra. Such a comparison is, 
however, extremely hampered by the fact that especially the soft 
X-ray emission is heavily absorbed by the interstellar medium. 
Another shortcoming of existing X-ray observations is their 
rather low spectral resolution. 

As an example, Fig.[T2]gives the computed X-ray flux from a 
nebular model as seen from a distance of 1 kpc, with and without 
interstellar absorption. The model used here is the one shown in 
the bottom panel of Fig. |4] It is of intermediate age, 5640 yr, 
and with steflar parameters of M = 0.595 Mq, L ^ 5200 Lq, 
and Teff - 7 1 500 K it is a typical representative of most objects 
with known X-ray emission. The X-ray spectrum is not expected 
to change much while the star evolves across the HR diagram 
since the characteristic temperature of the X-ray emitting region 
remains rather constant during the high-luminosity part of evo- 
lution (cf. Fig.[T8]in Sect.Ell). 

The upper panel of Fig. [T2]demonstrates the effect of (inter- 
stellar) extinction on the fully resolved X-ray spectrum, assum- 
ing a column density of 8x 10^" hydrogen atoms per cm^. The ex- 
tinction cross section per hydrogen at om as a function of X-ray 
energy was calculated according to I Morrison & McCammonI 
( I1983L Table 2). In the lower panel we show a simulation of the 
expected count rates as they would be measured by the XMM- 
Newton satellite with its limited spectral resolution, using an ap- 
propriate response matrijy describing the complex response of 
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Fig. 12. Top: Spectral flux densities computed from the model 
shown in the bottom panel of Fig. |4] (sequence No. 6a-HC2) 
at a distance of 1 kpc, non-attenuated (black) and attenuated 
by an intervening hydrogen column density of A^e = 8 x 



10 cm" (grey). Bottom: Simulated count rates per spectral bin 
of 0.03 keV width for the EPIC camera of the XMM-Newton 
satellite, with (gray) and without (black) absorption. 
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Fig. 13. Top: Spectral flux densities computed from a model se- 
lected from sequence No. 10a-HC2 at a distance of 1 kpc, non- 
attenuated (black) and attenuated by an intervening hydrogen 
column density of A^h = 6 x 10^' cm"^ (grey). The stellar param- 
eters are: M = 0.696 M©, L = 8403 Lq, T^sf = 199 064 K, t = 
695 yr. Bottom: Corresponding simulated count rates per spec- 
tral bin of 0.03 keV width for the EPIC camera of the XMM- 
Newton satellite, with (gray) and without (black) absorption. 



the energy channels of the EPIC camera to irradiation by X-ray 
photons of given energy. For comparison with published X-ray 
spectra of PNe, we have re-binned the XMM energy channels 
such that the full energy resolution of 5 eV is reduced to 30 eV. 
At this resolution, the only prominent feature seen at w 0.55 keV 
(22 A) belongs to the strong O vii complex. Other line features 
cannot be seen. 

A second example is shown in Fig.[T3]where the X-ray emis- 
sion was computed for a model from sequence No. 10a-HC2 
with M = 0.696 M©, and L ^ 8400 L©, T^s * 200 000 K, repre- 
sentative of, e.g., NGC 7027. There is more flux at high energies 
as compared to the case in Fig.[T2]because the bubble gas is now 
hotter (cf. Fig.|5] t op). 

According to iKastner et al.l (1200 ih . the absorption towards 
NGC 7027 is very high, A^h ~ 6 x 10^' cm"^, shifting the max- 
imum of the observed X-ray emission to higher energies (bot- 
tom panel of Fig. \T3[ . Our simulated XMM-Newton spectrum 
is very similar to that of NGC 7027 observed by the Chandra 
X-Ray Observatory (cf. lKastner et alj|200 1): the maximum flux 
occurs between 0.8 and 0.9 keV, probably due to the Ne ix com- 
plex. Indeed, by looking at the top panel of Fig. [T3] one sees a 
strong line blend around 0.9 keV (Ne ix) just at the maximum of 
the reddened flux density distribution. Ne ix is also strong for the 
model shown in Fig. [121 but far away from the maximum of the 
(reddened) flux density distribution. Because of the larger bubble 
temperature we see also peaks generated by Mgxi (^1.3 keV) 



and Si XIII (^1.9 keV). The Ovii complex which is dominant at 
lower temperatures (cf. Fig. \T2[ . is not visible anymore at this 
spectral resolution 0. 

We see from these examples that the flux detectable from 
PNe is quite limited to a certain energy range: at high energies, 
the flux decreases due to the physical structure of the bubble, and 
at low energies the emission is heavily absorbed by the interven- 
ing interstellar matter, and possibly also by dust and neutral gas 
in and around the PN itself. The range useful for detecting X- 
rays appears to be between ^ 0.2 and 2 keV (corresponding to 
the wavelength range 6 ... 60 A). 

5.1. Absorbed X-ray surface brightness distributions 

We have akeady discussed the intensity distributions of the X- 
rays emitted from PNe bubbles which are subjected to different 
treatments of the heat conduction in Sect. 14.21 Like the spec- 
tral energy distribution, also the radial intensity distribution is 
expected to depend on absorption. The reason is the large vari- 
ation of the absorption with frequency (see Figs. [12] and \T3[ . 
The emission from the cooler, denser bubble gas is more heavily 
absorbed than that from the hotter but less dense gas further in- 



' Note that all the computations presented in this work were per- 
formed with the same chemical abundances, viz. those listed in Table[T] 



1'10' 



S'lO'-^-^ 



6'10'-^- 



4'10''^- 



2'10''^- 



0.0 









, 






Intensity 5 -100 A 


NH 8 10'° cm" 








- 


■ 


1 , . . 



0.5 1.0 1.5 

Impact parameter p [10"cm] 



2.0 



1-10 



77 8-10-- 



6'10'- 



4'10'" 



2- 10"' 




No extinction 
NH= 8 10'° cm" 



intensity 5 - 1 00 A 



0.0 



0.5 




Impact parameter p [10"cm] 



1.5 



2.0 



Fig. 14. Top: Intensity profiles from Fig. [TO] middle panel, sub- 
ject to an absorption corresponding to a hydrogen column den- 
sity of A^H = 8 X 10^° cm"^ (grey)- Bottom: The same for the 
bottom panel of Fig. [TO] The black lines represent in both panels 
the unattenuated intensities. All intensities are integrated from 
5-100 A, corresponding to 2.5-0.125 keV. 



wards. Consequently, the emission is weighted towards the hot- 
ter, less dense inner regions of the bubble. 

This effect is illustrated in Fig. [14]where we show examples 
of radial X-ray intensity profiles with and without extinction for 
our two treatments of heat conduction. In both cases the centre- 
to-limb variation is considerably reduced if interstellar extinc- 
tion is applied, and the maximum emission is shifted inwards. 
With method 2 the X-ray intensity appears to fill the region in- 
side the nebular rim nearly homogeneously. 

The radial distribution of the X-ray intensity offers obviously 
a possibility to put constraints on the physics of thermal conduc- 
tion in very rarefied plasmas. Although it is difficult to estimate a 
definitive amount of center-to-limb variations from the existing 
X-ray maps, we note that the recent X MM-Newton observations 
of NGC 3242 presented by iRuiz et al.. (.2006.) show a smooth 
distribution of X-ray emission with no apparent limb brighten- 
ing. Thus we are inclined to conclude that our method 2 (labeled 
HC2) appears to be more appropriate to describe heat conduc- 
tion in PNe. We will therefore discuss in the following sections 
only this case. 

Although limb brightening is not really detectable, it is 
obvious that the X-ray intensity is often not homogeneously 
distributed across the bubble. In the heat conduction model 
this would be an indication that the conduction efficiency is 
not the same in all directions, maybe due to the presence of 
small-scale magnetic fields. It appears more likely, however, 
that the brightness variations are caused by non-uniform (intra- 
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Fig. 15. X-ray luminosity of the bubble for the energy interval 
0.45-2.5 keV (5-28 A) as a function of the stellar effective tem- 
perature as predicted by three model sequences with central star 
masses 0.565 Mq, 0.595 M©, and 0.696 M© (HC2 models only, 
see inset). The star symbols indicate the positions of the ob- 
jects from Table [3] The evolution times covered by the model 
sequences are the same as in Fig.fTTl 



ne bular/interst ellar) extinction across the nebula (see discussion 
in i Kastner et al. 2002). as found for many PNe in the visu al 
wavelength range (e.g. lSandin et al.ll2008HTsamis et al.ll2008l) . 



5.2. X-ray luminosities 

Now we compare global X-ray properties, in particular the lu- 
minosities, predicted by our models with the recent observa- 
tions made by the Chandra and XMM-Newton satellites. We dis- 
carded objects with a WC central star because wind and hot bub- 
ble have a hydrogen-poor composition which is completely dif- 
ferent from the (hydrogen-rich) composition assumed in our hy- 
drodynamical simulations. The X-ray properties of the remain- 
ing five objects are hsted in Table[3]and are taken from the com- 
pilation of . Guerrero et al.l (12005 a) in order to have a homoge- 
neous data set to compare with. Note that now the X-ray lumi- 
nosities refer to the energy interval of 0.45-2.5 keV (5-28 A) 
only in order to avoid large and uncertain corrections due to in- 
terstellar absorption at lower energies. The low energy cut-off 
used here leads to X-ray luminosities that are, in some cases, 
considerably lower than those quoted in the discovery papers. 
Additional changes are due to the adjusted distances used in our 
work. The X-ray data from the PNe are supplemented by the 
corresponding data of the stellar winds if available. 

The empirical data of Table [3] allow already some interest- 
ing conclusions: Compared with the stellar luminosity, the X-ray 
contribution from the shocked stellar wind (i.e. from the bub- 
ble), for the energy band considered here, is very small and is 
only about 10"^ i-star- Also the fraction of the stellar wind power 
that is converted into X-ray emission is quite small: only 1% to 
0.01%. Note that these values depend on the definition of the 
energy band, notably on the boundary of the low-energy region 
where most of the X-rays are emitted (see Fig.|9]l. 

We present in the following figures a detailed comparison 
of the objects listed in Table [3] with the predictions of our hy- 
drodynamical models with heat conduction treated according to 
method 2. We begin with the "X-ray Hertzsprung-Russell dia- 
gram" of Fig. [15] where we present the general evolution of the 
bubble's X-ray luminosity as a function of the stellar effective 
temperature for sequences with different central stars. The in- 



Table 3. Relevant parameters of objects with detected X-ray emission. The luminosity, Lx, in the wavelength band 5-28 A (0.45- 
2.5 keV), corrected for extinction and adjusted according to the distances used h ere, and typical values of temperature, Tx, and pres- 
sure, Px, of th e X-ray emitting volu me as determined bv lGuerreroet al.l (12005 al) . The effective temperatures are either spectroscop- 
ically derived dMendez et alJll992h or Hen Zanstra temper atures (NGC 6543 and NGC 7027: Gorny, priv. comm.). The distances 
are e ither spectroscopic one s (NGC 2392 and NGC 3 242: 'Pauld rach et al.ll20()4 ) or based on expansion parallaxes (NGC 6543: 
iReed et al. 1999: NGC 7027: ISchonberneret al.ll2005a ). T he distance of NGC 7009 is again from Gorn y (priv. comm.). Refe rences 
for the central-star wind data: (1) Pauldrach etal. 20041: r2) lTinkler & Lamersll2002l Table 1 therein; (3') lGeorgiev et alJl2008i Table 
1 therein. All objects listed in this table have central stars with normal, i.e. hydrogen-rich, composition. 



Object 



Teft- Distance 
(K) (kpc) 



(Lo) 



logLx 
(ergs"') 



Tx 
(K) 



^x 

(dyne cm" 



l0g(Lx/Asta,) 



log M„i„d 
(Moyr-') 



'^wind 

(kms"') 



^wind 

(L©) 



log(Lx/L„i„d) Ref. 



NGC 2392 


40000 


1.67 


5010 


31.20 


2.0 X 10'^ 


2.6 X lO'** 


-6.09 


-7.74 


420 


0.26 


-1.81 


(1) 


NGC 3242 


76000 


1.10 


3120 


30.90 


2.2 X 10* 


7.6 X 10-' 


-6.18 


-8.40 


2400 


1.89 


-2.96 


(1) 


NGC 6543 


67 000 


1.00 


1590 


31.00 


1.7x10" 


1.8x10-'* 


-5.78 


-7.73 


1340 


2.75 


-3.01 


(3) 


NGC 7009 


81000 


1.50 


3600 


31.34 


1.8 X lO*" 


1.4x10-'* 


-5.80 


-8.55 


2770 


1.78 


-2.49 


(2) 


NGC 7027 


200 000 


0.80 


6250 


31.12 


7.9 X 10* 


1.3x10-' 


-6.26 


— 


— 


— 


— 





crease of the X-ray luminosities with time seen in Fig.fTTItrans- 
lates into a corresponding increase with effective temperature. 

Given the distance uncertainty to individual objects, the 
agreement of our models with the existing observations is very 
good. The X-ray luminosity of all five objects from Table[3]can 
be explained by sequences with central stars between 0.6 and 0.7 
M0! We repeat that the run of the X-ray emission with time (or 
stellar effective temperature) depends on the stellar wind prop- 
erties and the expansion rate of the bubble, which also depends 
indirectly on the wind. The evolution of the X-ray luminosity as 
seen in Fig.[T5]is the result of the wind model shown in Fig.|2] 
(0.595 Mq). The sample of PNe with confirmed soft X-ray emis- 
sion is, however, too small to prove or disprove the wind model 
used in our simulations. 

Our models without thermal conduction fail completely in 
explaining the observations because their X-ray luminosities are 
too low by about two orders-of-magnitude at the positions of the 
observed objects in Fig.[T5](cf. also Fig.fTTTi. 

In order to avoid possible systematic errors caused by dis- 
tance uncertainties we relate in Fig. [16] the X-ray luminosities 
to both the stellar luminosity, Lx/i^stai, and wind luminosity, 
Lx/i'wind- The basic result found from Fig. [15] is confirmed in 
Fig.[T6](top panel): our models with thermal conduction accord- 
ing to method 2 (HC2 models) predict soft X-ray emission in 
good agreement with the observations (except for NGC 7027, 
see discussion below). There appears to be a slight preference 
for our 0.7 M0 models, but considering the uncertainties of Lx 
we don't think that this is a real effect. 

The lower panel of Fig. [16] corresponds to that shown in the 
bottom panel of Fig. [TT] although the difference between the 
tracks appears somewhat larger The reason is the smaller en- 
ergy range considered: the temperature sensitive emission at the 
high energy end has now more weight. Again we have a satisfac- 
tory agreement between theory and observations. Less than 1 % 
of the wind power is converted into X-ray emission. Note, how- 
ever, that the mass-loss rates which enter into the wind power 
are notoriously uncertain and may vary by up to a factor of ten 
between different authors! 

As mentioned above, NGC 2392 has an unusually slow wind 
for its position in the Hertzsprung-Russell diagram compared 
with the other objects listed in Table[3] giving it the lowest wind 
luminosity and thus the highest LxlL^md ratio of all objects from 
the sample. Its X-ray luminosity agrees with our models (cf. Fig. 
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Fig. 16. Top: X-ray luminosity of the bubble (0.45-2.5 keV) 
over stellar luminosity vs. effective temperature. The model se- 
quences are the same as in Fig. [15] (see inset). Bottom: X-ray 
luminosity of the bubble (0.45-2.5 keV) over stellar wind power 
vs. stellar effective temperature. The observations from Table [3] 
are again plotted as stars. Note that NGC 7027 cannot be shown 
in the bottom panel because its stellar wind properties are un- 
known. 



[T6] top), but its wind power is too small, and hence Lx/L^md 
too large (cf. Fig. [16] bottom). It appears possible that fast out- 
flows/jets co ntribute to the X-ray e mission of NGC 2392 (see 
discussion in lGuerrero et al.ll2005bl) . 
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Fig. 17. Lx/Lstai vs. wind speed, again for our three model se- 
quences from Fig. [15] and the objects from Table [3] (except for 
NGC 7027 for which the central star wind is not known). The 
X-ray range is again 0.45-2.5 keV. 

A further interesting comparison is shown in Fig. [TT] where 
Lx/Lscdr is plotted over ywind- This figure demonstrates clearly 
that our models, thanks to the inclusion of heat conduction, give 
a fully consistent description of the observed X-ray luminosities 
also in terms of the observed large wind speeds! Although the 
wind of NGC 7027 is not known, it must exist because we ob- 
serve the X-ray emission from the shocked wind gas. Judging 
from the central star's position close to the white dwarf domain, 
we estimate a wind speed close to the right boundary of Fig. [17] 
viz. of about 7000. . . 8000 km s"'. 

5.3. X-ray temperatures 

A stringent test of our models is the determination of a charac- 
teristic X-ray emission temperature, Tx, which can be compared 
with the measurements. We computed Tx from our models by 
weighting the electron temperature Tg{r) within the bubble (ri- 
ra) with the volume emissivity integrated over the respective en- 
ergy range, E1-E2 (0.45-2.5 keV or 5-28 A): 



47r r 

^X Jr, 



r'T,(r)T]x(r)dr, 



where Lx is the X-ray luminosity. 



Lx=4 



and 






r T]x(r)dr, 



-£2 



'7x(r)= r r](T,(r),n,(r),E)dE, 
Jei 



(17) 



(18) 



(19) 



is the volume emissivity in the energy band E1-E2. 

The result is seen in Fig.[T8](top panel), covering the whole 
evolution from the onset of the bubble formation until the white 
dwarf stage is reached (cf. Fig. [TTT i. In general, Tx increases 
rapidly to above 10^ K while a hot bubble is formed beyond the 
wind shock. Then heat conduction becomes effective, and Tx in- 
creases only slowly or remains nearly constant until maximum 
wind power is reached. This refers mainly to the 'evaporation' 
phase during which the the conduction front advances outwards 
(see Sect. [3]). Afterwards, Tx drops in line with the wind power 
to about 10*' K ('condensation' phase). 

The maximum X-ray temperature achieved during the evo- 
lution across the HR diagram depends on the central star mass: 
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Fig. 18. Top: LxlL^tm for the 0.45-2.5 keV range as a function 
of the characteristic temperature of the X-ray emitting region, 
Tx, again for the three sequences from Fig. [15] The individual 
'tracks' cover the same age spans as shown in Fig.[TT] Note that 
the observed values of Tx (from Table[3]l have been derived from 
the best fit to the spectral energy distributions, while the Tx for 
a given hydrodynamical model have been computed according 
to Eq. ( [TTI i. Bottom: Tx as a function of the wind speed, Wwind, 
computed for the three model sequences from Fig. [15] and Tq 
according to Eq. ([T2l l (dashed). 



Thanks to its powerful wind, the sequence with the 0.696 M0 
model reaches a much larger Tx than the sequence with the 
0.565 Mq star, viz. 4.5 x lO*- K instead of only 1.6 x 10^ K. 

Assuming that our computed characteristic X-ray tempera- 
tures are reasonable representatives of the ones derived from ob- 
servations, we find excellent agreement between theory and ob- 
servations, except for NGC 7027 whose Tx (8 x 10^ K) is much 
larger than our models predict (see also below). The remaining 
objects are rather well matched by the 0.595 Mq models with 
heat conduction (sequence No. 6a-HC2). 

The bottom panel of Fig. [T8]illustrates how the temperature 
of the X-ray emitting region is modified by thermal conduction 
from a value close to the one predicted by Eq. ([T2I 1. At very low 
wind speeds, typical for the 'early-wind' phase, thermal conduc- 
tion is unimportant, and Tx equals the post-shock temperature 
(dashed line in the bottom panel of Fig. [TSl l. At faster stellar 
winds (and larger post-shock temperatures), Tx levels off as a 
consequence of heat conduction and becomes fairly independent 
of the wind speed. 



5.4. Correlations with bubbie radius 

For a given central star mass, our models predict an increase of 
the X-ray luminosity with time (bubble radius) during the main 
phase of evolution (see Sect. l4.3l Fig.[TTl top). On the other hand, 
the same models predict a decrease of Lx with bubble radius 
at given central star effective temperature. This is because, for 
given reft, Lx increases with stellar mass (see Sect. 15. 21 Fig.fTSll. 
while the bubble size decreases due to the shorter evolution time 
scales of the more massive stars. Likewise, in the presence of 
heat conduction the mean temperature of the hot bubble (~ Tx) 
depends on both the wind power and the bubble radius according 
to Eq. (fT4l i. but the dependence is weak. Hence, the expected cor- 
relations of Lx and Tx with bubble radius 7?hb depend critically 
on the mass distribution and ages of the observed objects. 

In this contex t we have to discuss the recent findings by 
iKastner et al.l (2008) that X-ray temperature and luminosity ap- 
pear to decrease with bubble radius. Apart from the fact that 
these results suffer from uncertainties of the distances of the in- 
dividual objects, which the authors do not consider at all, we 
point out that: 

1 . Figure 4 of iKastner et al.l (l2008l) contains a mix of objects: 4 
PNe have hydrogen-deficient and 3 PNe have hydrogen-rich 
central stars. The objects from both groups have certainly /jer 
se distinct physical properties and also different evolutionary 
histories. 

2. The 3 objects with a hydrogen-rich central star (NGC 2392, 
NGC 6543, NGC 7009) have about the same bubble radius 
of~0.1pc. 

The claimed anti-correlations with bubble radius are thus proba- 
bly purely artificial, based entirely on an inappropriate combina- 
tion of objects with different evolutionary background. A much 
larger, homogeneous sample of objects is certainly necessary to 
construct trustworthy correlations between observable quantities 
that can be compared with theoretical predictions. 

5.5. Individual objects 

NGC 3242. This PN is well suited for a closer comparison 
with our models because its bubble is most likely spherical, as 
judged from the ring-like appearance of the rim. The position of 
NGC 3242 is very close to the 0.595 Mq track in all the previ- 
ous Figures ( fTsl [T6l [TTl andfrsTl . This implies that the evolution of 
wind power and X-ray emission as predicted by our model sim- 
ulation with heat conduction according to method 2 reflects the 
real situation in NGC 3242 surprisingly well. A detailed com- 
parison between the observed parameters of NGC 3242 and two 
models taken from sequence No. 6a-HC2 which embrace the 
observed position of NGC 3242 in the figures is presented in 
Table |4] The density and temperature structure of model 1 is 
shown in Fig. |4] (bottom panel). 

Given the uncertainty of the observed data and the fact that 
we did not attempt to make any fits to the observations, the 
agreement is very good, especially for the X-ray related data. 
Notice the nearly equal pressures and temperatures of the X-ray 
emitting gasQ Also mass-loss rates and wind speeds are, within 
the known uncertainties, in reasonable agreement. 

Quantities that depend directly on distance or distance 
squared differ by larger amounts. Our values for NGC 3242 
quoted in Table |4] are based on a distance of 1 . 1 kpc, accord- 
ing toUauldrach et al. (2004). Kudritzki et al. (2006) arrived at 

* The hot bubble is isobaric despite of the radial temperature and 
density gradients (cf. Fig.|4j. 



Table 4. Relevant parameters of two nebular models along the 
0.595 M0 track (sequence No. 6a-HC2) embracing the position 
of NGC 3242 in the HR diagram, compared with the observed 
properties of NGC 3242. For both observation and models, the 
X-ray data refer to the 0.45-2.5 keV energy band. 





Model 1 


Model 2 


NGC 3242 


Ref. 


M/Mq 


0.595 


0.595 


0.53, 0.63 


(1), (2) 


AstaiV-^O 


5 205 


5 051 


3162 


(1) 


T,« (K) 


71667 


80457 


75 000 


(1) 


?post-agb (yr) 


5 642 


6121 


^2 800 


(3) 


M„ind (Moyr"') 


9.7 X 10-*^ 


8.6 X 10-" 


4 x 10-" 


(1) 


!)„i„d(kms-i) 


2115 


2490 


2400 


(1) 


^x/^star 


5.0 X 10"^ 


6.8 X 10-' 


6.6 X 10-' 


Tabled 


^x/^witid 


7.2 X 10-* 


7.7 X 10-" 


1.1 X 10-^' 


Tabled 


TxCK) 


2.1 X lO" 


2.1 xlO" 


2.2 X 10" 


(4) 


Px(dyncm-^) 


9.0 X 10-'' 


8.0 X 10-" 


7.6 X 10-" 


(4) 


Prim (dyncm--) 


^8x 10-' 


=^7x10-" 


8.6 X 10-" 


(4) 


«rira (cm"^') 


^2900 


::=2400 


2 600 


(4) 


Oshcu (cm-^) 


^1200 


^800 


800 


(4) 


!>,im(kms-') 


14.0 


14.5 


19.5 


(5) 


Usheii (kms-') 


27.0 


28.0 


35.7 


(5) 



(l)IPauldrach et al."2004; 

(2) Kudritzki et al. 2006; 

(3) Kinematic age from Corradi et 
scaled down to the distance of 

(4) Ruiz et al. 2006; 
(5) .Schonbemer et ani2005bl . 



al.'2003', 
1 . 1 kpc used 



here; 



a larger distance, li = 1 .8 kpc, hence luminosity and kinematic 
age are larger: L - 7760 Lq with an age of about 4 600 yearfl 
These authors give also a larger mass-loss rate, 8 x 10"'' Mq yr"' , 
from the central star. 

Taken at face values, the numbers fisted in Table|4] indicate 
that the thermal pressure of the rim, Prim, exceeds that of the bub- 
ble by a small margin. In the model we have just the opposite sit- 
uation. The case of NGC 3242 could thus indicate that the wind 
power of the central star has achieved its maximum value al- 
ready below the present effective temperature of 75 000 K. This 
would be in line with the more recent wind computations con- 
ducted by lPauldrach et al.l (l2004l) according to which central-star 
mass-loss rates and wind power reach a maximum around stel- 
lar temperatures of 50 000 K, after which the thermal pressure 
of t he bubble drops below th at of the rim (for more details see 
Stef fen & Schonberneij|2006l Fig. 6). In our present simulations 
which are based on the older lPauldrach et al.l (11988) recommen- 
dations the wind power peaks later, close to maximum effective 
temperature (cf. Fig.|2l bottom left). 

One may then ask whether also other structures of NGC 3242 
are adequately represented by our models. First of all, the model 
has a dense rim and a large but less dense (attached) shell, and 
size and density ratios between rim and shell are about two and 
three, respectively. These values compare favorably with the ob- 
servations, and also the (electron) densities are in good agree- 
ment: Ruiz et al. (2006) derived electron densities for NGC 3242 



' Kinematic ages based on physical size and expansion velocity 
are problematic since the y depend strongly on the method used (see 
ISchonbemer et al.]l2005ah . The relatively large shell expansion speed 
indicates a significant acceleration during the previous evolution which 
lead to an underestimation of the age. 
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Fig. 19. To/j: normalised intensity distributions of NGC 3242 in 
[O m] and He ii from HST monochromatic images (F502N and 
F469N, respectively). The cuts are taken along the minor axes 
and scaled to the model sizes. The spikes at r = are due to 
the central star. Bottom: normalised intensity distributions of an 
appropriate hydrodynamical model with a 0.595 Mq central star 
that matches closely the observations for NGC 3242 (model 2 of 
Tab.H. 



of 2600 cm ^ in the rim and 800 cm ^ in the shell, respectively 
(cf. TableH. 

Figure[T9]gives a further illustration of the close relationship 
between NGC 3242 and our models: it shows the monochro- 
matic brightness distributions in two important emission lines 
for NGC 3242 and for an appropriate model very close to the 
observed HRD position of NGC 3242. It is remarkable that also 
the ionization structure in both the real PN and the model is such 
that only the rim is doubly ionized in helium. Yet the model is 
not perfect: According to Table|4]the (spectroscopic) expansion 
velocities of rim and shell are lower than the observed ones. The 
velocity differences, however, are well matched. 

NGC 7027. We note from Fig. [16] (top panel) that, although 
the X-ray luminosity of NGC 7027 compares well with those of 
the other objects shown in this figure, it is about 1 .5 dex below 
the model prediction. Also, the X-ray emitting region is con- 
siderably hotter, viz. with 8 x 10^ K about twice as hot as the 
maximum reached by our 0.696 M0 models of sequence No. 
10-HC2 (see Fig.[T8ll. We have verified that NGC 7027 is em- 
braced by our sequences No. 10 (no conduction) and No. 10a- 
HC (method 1), so we conclude that thermal conduction may 
still work but at a lower level than our method 1 predicts. A 
possible solution would be the presence of a weak magnetic 
field which suppresse s thermal conduction to some extent (cf. 

lBorkowskiet"ai]ll990l). 

Recently, ISabin et aP (l2007h reported indeed the detection 
of polarization across NGC 7027 by means of SCUBA obser- 
vations. The orientation of polarization indicates the presence 
of a toroidal magnetic field along the equatorial plane. Across 
the central cavity, no clear polarization is visible, thus any state- 
ments about the orientation of a possible magnetic field inside 
the bubble is impossible. If this interpretation for the rather low 
X-ray luminosity is correct, we must infer that heat conduction, 
in the particular case of NGC 7027, is not fully suppressed by 
the presence of a (weak) magnetic field. The field geometry is 



likely to play a role. Indeed, it see ms that the X-ray em ission is 
suppressed in the equatorial plane (iKastner et al.ll2002l) . 

Alternatively, the low X-ray luminosity of NGC 7027 might 
be related to t he pre sence of bipolar collimated outflows found 
bv'Cox e t alj (|2002|) . It is conceivable that a significant amount 
of potentially X-ray emitting matter is lost through multiple 
openings in the skin of the hot bubble created by the jets. 

NGC 2392. This is a rather peculiar object which poses a prob- 
lem for our models, as already noted above. Its central star has a 
wind speed that is much too low for the object's position in the 
Hertzsprung-Russell diagram at about 40000 K effective tem- 
perature. The wind speed is more typical for the end of the early- 
wind phase, although the mass-loss rate appears to be rather nor- 
mal (see TableO. This low wind speed is responsible for a wind 
luminosity which is about a factor of ten below the wind powers 
of the other objects of this study (cf. Fig. [3]). 

Despite of this, the X-ray emission from the hot bubble of 
NGC 2392 compares well with that of the other objects (cf. Fig. 
[T6] l and is much too high for an early-wind with speeds below 
500 km s"', as is predicted by our models (Fig.fTTb. It is possible 
th at additional X-ray e mission is provided by jets as proposed 
bv lAkashi et"an (l2008h . 

5.6. The UV emission lines 

Due to the steep temperature gradient across the bubble/nebula 
interface, there exists only a very narrow region which is suitable 
for the emission of UV lines from highly ionized species. We re- 
call here that our code computes the ionization of all elements 
considered (see Table [T]) time-dependently for the whole com- 
putational domain, i.e. also in the freely streaming and shocked 
stellar wind. As a byproduct of our simulation, we are thus able 
to compute al so the line emi s sion fr om the bubble/PN interface. 

Recently, iGruendl et alJ (l2004l) reported the detection of 
Ovi /1 1032 and 1038 A emission lines in FUSE spectra of 
NGC 6543. Since the central star is not very hot, the ioniza- 
tion within the nebular shell is too low as to account for O^^. 
The authors concluded that the O vi lines must originate from 
the conduction front at temperatures of about a few times 10^ K. 

Test calculations showed that the O^^ shell or 'pocket' is ex- 
tremely thin, only about 1 x 10'^ cm thick, which is comparable 
with the spatial resolution of our numerical mesh. We thus re- 
computed the sequences Nos. 6a, 6a-HC, and 6a-HC2 with a 
finer mesh (Ar = 3 x 10''* cm) for r < 1.5 x 10'^ cm in order to 
achieve a better resolution of the O^^ layers. 

Figure |20] illustrates how the ionization of oxygen varies 
within the hot bubble, and how the distribution of O^^ depends 
on the physical treatment of this region. Note that, since the 
thermal structure of the conduction front/contact discontinuity 
does not change much during the lifetime of the PN, the ion- 
ization fractions remain rather stationary with respect to the 
front/discontinuity. We find that the thermal structure close to 
the conduction front where O^^ prevails does not differ much 
between our two treatments of heat conduction. 

The physical conditions change rapidly across the conduc- 
tion front: ahead photoionization by the stellar radiation field is 
the dominant heating and ionization mechanism. For the stellar 
temperature shown in Fig. |20l O^^ and O^^^ are practically the 
only representatives of oxygen. Behind the front the tempera- 
tures are so high (<: 10^ K) that electron collisions determine the 
ionization state of the gas: now we have a mixture of O*"^, O^^, 
and O**^ in proportions ruled by the electron temperature. For 
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Fig. 20. Radial profiles of the ionization fractions of oxygen 
in the vicinity of the contact discontinuity /conduction front for 
different physical treatments: Top: without thermal conduction; 
middle: with thermal conduction according to method 1 ; bottom: 
with thermal conduction according to method 2. Thick lines re- 
fer to the nebular region, thin ones to the bubble region behind 
the contact surface/conduction front. The very thin O^^ 'pocket' 
is shaded for clarity. Details about the three models are given in 
Table 121 



instance, in the case without heat conduction the bubble temper- 
ature is so large that O*^^ (the bare oxygen nucleus) is the main 
constituent already right behind the contact discontinuity. Heat 
conduction lowers the temperature gradient behind the front, and 
consequently O''^ and O''^ are the main constituents of oxygen 
throughout a large fraction of the bubble. 

The transition between the photo-heated nebular gas (T^ ^ 
10^ K) and the shock-heated wind gas (Tg > 10^ K) occurs 
very abruptly across the contact discontinuity if there is no ther- 
mal conduction, and consequently only little room is left for the 
existence of O^^, simply because the temperature is either too 



low (nebula) or too too hot (bubble). With thermal conduction 
included, the temperature increases somewhat more gently and 
thus allows for a larger amount of O^^. This is also reflected in 
the total emission of the O vi lines listed in Tab.|5]for the models 
shown in Fig.|20] We point out that the computed O vi luminosi- 
ties are strongly fluctuating from model to model due to the poor 
numerical resolution of the conduction front. The O vi line fluxes 
listed in Tab.|5] are therefore appropriate averages over several 
contiguous models. 

According to the results for the heat conduction models in 
Tab.|5l the total flux in the Ovi Hne at A 1032 A emitted by the 
O^"^ pocket, assuming a distance of 1 kpc, would be roughly 
4 X 10"''' ergcm^^s"' . This is consistent with the estimates of 
Gruendl et al. (2004) for NGC 6543, although their FUSE mea- 
surements are based on a slit across the central cavity. The total 
O VI /i 1032 A line flux is certainly larger than the quoted value 
of ^2 X 10"'^ ergcm^^s"'. 

We conclude that the luminosity in UV O vi lines is gener- 
ated in a very thin transition layer which is not very sensitive to 
the effects of thermal conduction. In contrast, the X-ray luminos- 
ity depends sensitively on the efficiency of thermal conduction. 
In our models with heat conduction included, the O vi luminosity 
is comparable to the X-ray emission coming from more extended 
parts of the bubble (cf. Tables[3]and|5]l. 



6. Summary and conclusions 

We presented a detailed numerical approach towards an under- 
standing of the diffuse soft X-ray emission from planetary neb- 
ulae based on the concept of thermal conduction. Since thermal 
conduction is a physical process inherent to all hydrodynami- 
cal systems, becoming important wherever the mean free path 
of the electrons is sufficiently large, we included a thermal con- 
duction module into our ID radiation-hydrodynamics code. We 
were able to compute the thermal structure of the shocked wind 
gas inside the nebular cavity self-consistently with the hydro- 
dynamics once the stellar AGB remnant, its initial circumstellar 
envelope, and the post- AGB wind model are specified. By doing 
so, there is some freedom in the treatment of heat conduction in 
cases where the mean free path of the electrons becomes com- 
parable to the characteristic temperature scale length. Magnetic 
fields that may play a role for the shaping of PNe are not consid- 
ered. 

Thermal conduction has a twofold effect favoring thermal X- 
ray emission from the shocked wind gas: it lowers the temper- 
ature gradient across the bubble/nebula interface and heats the 
cool nebular matter, forcing it to 'evaporate' inwards. With time, 
heat conduction accumulates additional matter in the bubble 
with a characteristic temperature of some 10* K, which quickly 
dominates the bubble's mass budget. The amount of added mat- 
ter is controlled by heat transfer across the bubble from the in- 
verse wind shock towards the nebula. The bubble mass may in- 
crease or decrease with time, depending on the evolution of the 
stellar wind power. 

We s elected several of the hydrodynamic sequences pre- 
sented in lPaper j and recomputed them, without changing any of 
the other parameters or boundary conditions, with our heat con- 
duction treatment included. The X-ray emission was computed 
post-facto by employing the CHIANTI code, slightly adapted 
to match our purposes. We were able to find good agreement 
with existing observations, both with respect to X-ray luminos- 
ity and surface brightness distribution, if method 2 is used for the 
treatment of heat conduction. At the same time, the stellar wind 



Table 5. Dependence of the luminosity of the O vi lines at A 1032 A and 1038 A on the treatment of heat conduction for models 
similar to those shown in Fig.|4] but computed with higher spatial resolution. L(li/3) and L(Oiii 5007 A) refer to the H/3 and [Oiii] 
luminosities of the whole nebula. 



No. 


M 
(Mo) 


L 

(Le) 


T-eff 

(K) 


Thermal 
conduction 


(Le) 


L(0 


III 5007 A) 

(Le) 


L(Ovil032A) 
(Le) 


L(Ovil038 A) 
(Lq) 


6a-HR 

6a-HCHR 

6a-HC2HR 


0.595 
0.595 
0.595 


5 260 
5 297 
5 292 


67 978 
65 160 
65 541 


no 

method 1 
method 2 


21.7 
23.9 
23.6 




168.8 
174.0 
173.3 


4.41 X 10"' 
11.3x10-3 
10.1x10-3 


2.21 X 10-3 
5.63 X 10-3 
5.04 X 10-3 



model employed here, based on the theory of radiation-driven 
winds, implies wind speeds in excess of 1000 kms-' for typi- 
cal nebular models (cf. Fig.fTTli. We conclude that our modeling 
is fully consistent with both the observed wind speeds of more 
than 1000 kms-' and the observational evidence that the soft 
X-ray emission comes from regions with electron temperatures 
of about 2 X 10^ K (see also Table [3]l. In contrast, our models 
without heat conduction fail to reproduce the characteristics of 
the observed X-ray emission by large amounts if realistic wind 
parameters are assumed. 

The basic findings from our modeling can be summarized as 
follows; 



- The energy budget of the hot bubble is substantially altered 
by radiative losses at the bubble/nebula interface, both in 
models with and without thermal conduction. 

- Heat conduction has a strong influence on the thermal struc- 
ture of the hot bubble, but leaves the dynamics of the whole 
system is virtually unchanged, i.e. the shaping of the nebular 
shells and their expansion properties are not affected. 

- In the absence of heat conduction, the temperature of the hot 
bubble depends only on the velocity of the central star wind. 
The resulting X-ray temperature is much too high, and the 
X-ray luminosity much too low, to explain the existing ob- 
servations. 

- If heat conduction is substantial, the bubble temperature is a 
function of the stellar wind power and the bubble size. The 
compact nebulae around massive, short-lived central stars 
are therefore expected to have hotter bubbles than the neb- 
ulae of low-mass, slowly evolving central stars. 

- The X-ray luminosity is determined by the total emission 
measure, which decreases by expansion and increases by 
'evaporation' from the main nebula. The competition be- 
tween expansion and 'evaporation' rules the temporal evo- 
lution of the X-ray luminosity. 

- According to our PN simulations with a time-dependent 
wind model, the X-ray luminosities increase with time dur- 
ing the main phase of evolution when the wind power in- 
creases with time, in contrast to the case of a constant wind 
power that would imply a decrease of the X-ray luminosity 
with time. 

- For the energy range usually considered for planetary neb- 
ulae, the X-ray luminosity is below about 1 % of the stellar 
wind power, and between 10-"* and 10-** of the stellar bolo- 
metric luminosity. The exact numbers depend on the actual 
evolutionary state and the wind model used. 

- For given effective temperature, the X-ray emission is largest 
for nebulae around massive central stars even though their 
small bubbles contain only little mass. The X-ray emission 
measure of their bubbles is nevertheless large, primarily be- 
cause the high electron densities overcompensate the small 



amount of X-ray emitting gas that massive central stars can 
accumulate during their short lifetimes. 

- In contrast to the X-ray luminosity, the emission of the UV 
O VI lines is not very sensitive to the effects of thermal con- 
duction. 

- Magnetic fields must be absent or extremely weak in all ob- 
jects with diffuse X-ray emission since their presence would 
strongly depress thermal conduction and hence 'evapora- 
tion' . The absence or weakness of magnetic fields implies 
also that they cannot be responsible for shaping these ob- 
jects. 

Our numerical treatment as described here leads to X- 
ray luminosities that are substa ntially below those found by 
IZhe kov & Perinottol(ll996lll998h with their analytical approach, 
although their wind model and central-star evolution are very 
similar. We believe that their analytical approach overestimates 
the evaporation rate because radiative cooling of the gas at the 
conduction front by line emission is not considered (see also 
Sect.O. 

It is likely that the bubble/PN interface becomes dynamically 
unstable, leading t o direct mixing betw een hot bubble and cool 
nebular matter (cf. lStute & Sahaill2006l) . The net effect would be 
similar to heat conduction, i.e. a reduction of temperature gradi- 
ents allo wing a more efficie nt X-rays emission, as already men- 
tioned by IChu et al] ( 1997|). The 20 simu lations of a spherical 
bubble performed by lStute & Sahail (l2006h suggest that the mix- 
ing region is confined to a rather thin shell at the surface of the 
bubble. If so, the X-ray emission would be limb brightened, very 
similar to our models computed according to method 1 . For the 
time being, the existing observations do not seem to be consis- 
tent wdfli_suchajimb brightening. 

iGeorgiev et al.l (1200 8) found recently that the wind of the 
central star of NGC 6543 is less deplete d in iron compared t o 
the plasma emitting the diffuse X-rays dGeorgiev et alJ 120061) . 
They concluded that the X-ray emitting plasma "is derived from 
nebular gas rather than the stellar wind". This finding is fully 
consisting with the heat-conduction models. Clearly, more stud- 
ies of this kind would be very rewarding. 

Finally we want to emphasise that the models introduced 
in this work are entirely based on normal chemical composi- 
tions in the stellar and circumstellar envelopes. The results pre- 
sented here should therefore not be used to interprete X-ray 
emission from objects with hydrogen-deficient central stars such 
as BD-h30°3639, NGC 40 because their evolution is different 
and not yet understood. Additionally, one has to deal with heat 
conduction and X-ray emission in a practically hydrogen-free 
plasma. 
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